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ABSTRACT 


An  experimental  apparatus  was  designed,  constructed 
and  tested  for  the  measurement  of  very  low  partial  pressures 
of  H2S  and  CC>2  in  aqueous  ethanolamine  (EA)  solutions  at 
high  temperatures. 

The  apparatus  was  tested  by  measurements  of  CC>2  and 
H2S  solubility  in  the  region  where  previous  data  were  avail¬ 
able.  Tests  were  carried  out  for  C02  at  60°,  80°  and  100°C 
and  monoethanolamine  (MEA)  concentrations  of  2.5  and  5.0 
normal  (N) .  For  H2S  measurements  were  made  at  80°  and  100°C 
with  MEA  concentrations  of  2.5  and  5.0  N.  Good  agreement 
with  the  literature  proved  the  reliability  of  the  apparatus. 

Measurements  were  also  made  at  conditions  where  no 
previous  data  were  available.  The  80°  and  100°C  isotherms, 
for  both  pure  CC>2  and  H2S  in  2.5  N  MEA,  were  extended  to 
very  low  partial  pressure  of  the  acid  gas  (about  0.0001 
psia) .  Measurements  for  mixtures  of  H2S  and  C02  in  5N  MEA 
at  100°C  were  also  made  for  very  low  loadings  of  C02  and 
H2S  (each  being  less  than  about  0.1  mole/mole  of  MEA)  in 
the  liquid  phase. 

13 

The  theoretical  models  proposed  by  Klyamer  et  al. 
and  Kent  and  Eisenberg^  were  used  to  predict  partial  pres¬ 
sures  of  C02  and  H2S.  Modifications  were  carried  out  in 
the  model  of  Klyamer  et  al.  by  incorporating  fugacities  of 
CC>2  and  H2S.  The  functions,  for  activity  coefficients  of 

iv 
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ions  in  the  solution,  were  modified  to  include  the  effects 
of  temperature  and  initial  EA  concentrations.  Predictions 
from  these  three  models  (namely  Klyamer  et  al. ,  Kent  and 
Eisenberg,  and  modified  Klyamer  et  al.)  were  compared  with 
the  experimental  data  obtained  in  this  work,  as  well  as 
that  from  the  literature.  For  most  cases  the  Modified 
Klyamer  et  al.  model  gives  the  best  fit  of  the  experimental 
data.  For  the  simple  systems  (i.e.,  aqueous  EA  solutions 
containing  either  CC^  or  K^S) ,  however,  the  model  of  Kent 
and  Eisenberg  gives  a  slightly  better  fit  of  the  experi¬ 
mental  data. 
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NOMENCLATURE 


a  Ratio  between  the  un-ionized  ethanolamine 

activity  and  the  activity  of  water 
a  Constant  in  Redlich-Kwong  equation 

b  Constant  in  Redlich-Kwong  equation 

f  Fugacity 

H  Inverse  of  Henry's  Law  constant 

K  Equilibrium  constant 

k. .  Characteristic  constant  for  i-j  interaction 
il 

m  Molality  of  ethanolamine  in  Klyamer  et  al. 

correlation;  molarity  of  ethanolamine  in  Kent 

and  Eisenberg  correlation 

N  Normality  of  ethanolamine 

P  Critical  pressure 

o 

PC0  C02  Part^al  pressure 

PR  s  H2S  partial  pressure 

R  Gas  constant 

T  Temperature 

T  Critical  temperature 

c 

v  Molar  volume 

Vc  Critical  volume 

XCQ  Moles  of  C02/mole  of  ethanolamine 

2 

XTT  _  Moles  of  H^S.mole  of  ethanolamine 
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Mole  fraction  of  component  i 
z  Compressibility  factor;  carbamate  ion 

concentration 

zc  Critical  compressibility  factor 

Electrical  charge  on  an  ion 

Greek 

a  Activity  coefficient  of  water 

y  Activity  coefficient  of  an  ion 

y  Ionic  strength 

<J>  Fugacity  coefficient 

ft  Dimensionless  constant  in  Redlich-Kwong 

equation 

w  Acentric  factor 
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INTRODUCTION 


Aqueous  ethanolamines  are  widely  used  in  industry 
for  removing  f^S  and  C0£  from  gas  streams.  Vapor-liquid 
equilibrium  solubility  data  for  I^S  and  C02  and  their 
mixtures  in  EA  solutions  are  needed  to  aid  in  the  design  of 
amine-type  natural  gas  sweetening  units.  Data  of  this  type 
allow  the  estimation  of  the  equilibrium  capacity  of  the 
amine  solution  for  each  different  sweetening  system  design 
for  the  particular  conditions  of  temperature,  pressure,  and 
concentration  found  in  the  processing  units.  The  data  may 
be  used  to  check  the  feasibility  of  more  economical  designs 
not  yet  tested  by  plant  experience  and,  therefore,  provide 
a  means  to  guide  a  program  intended  to  reduce  sweetening 
plant  investment  and  operating  costs. 

Solubility  data  for  acid  gas-aqueous  EA  systems 
have  been  collected  by  numerous  investigators  over  the  last 
forty  years.  However,  the  data  published  to  date  do  not 
cover  all  the  necessary  ranges  of  solution  strengths  and 
temperatures.  The  region  of  low  loadings,  in  terms  of 
moles  of  acid  gas/mole  of  EA  and  low  acid  gas  partial  pres¬ 
sure  (where  most  commercial  units  operate) ,  has  not  been 
sufficiently  covered.  This  project  was  undertaken  to  cover 
this  important  region  of  very  low  partial  pressures  of  acid 
gases  and  high  temperatures.  To  make  these  measurements  a 
low  pressure  apparatus  had  to  be  designed  and  constructed. 
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To  be  able  to  use  this  apparatus  with  confidence  for  regions 
where  no  data  exist,  some  measurements  had  to  be  made  to 
reproduce  existing  data. 

While  the  equilibrium  solubility  data  are  useful  in 
themselves,  it  is  desirable  to  have  a  computer-oriented 
correlation  which  reproduces  the  experimental  results,  and 
which  can  be  extrapolated  outside  the  range  of  the  data 
with  some  degree  of  confidence.  So  far  only  two  solubility 
models  have  been  proposed  for  CC^-H^S-water-EA  systems  by 
Klyamer  et  al.^  and  Kent  and  Eisenberg^ .  The  starting 
point  for  both  of  these  models  is  the  postulation  of  chem¬ 
ical  reactions  occuring  in  the  aqueous  EA  solution  contain¬ 
ing  the  acid  gases.  The  model  of  Klyamer  et  al.  uses  activ¬ 
ity  coefficients,  for  ions  in  the  solution,  to  account  for 
the  non- ideality  of  the  liquid  phase.  Kent  and  Eisenberg 
have  incorporated  all  the  non-idealities  in  two  empirical 
parameters.  They  determined  these  parameters  from  experi¬ 
mental  data.  Neither  of  these  models  is  reliable  over  all 
conditions  of  solution  strengths  and  temperatures.  Modifi¬ 
cations  in  the  Klyamer  et  al.  model  (made  in  this  work) 
improved  substantially,  the  reliability  of  prediction  of 
equilibrium  partial  pressures  from  a  knowledge  of  liquid 
composition  and  temperature. 

In  Chapter  I,  design  considerations,  construction, 
and  operating  details  of  the  experimental  apparatus  are 
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described.  In  Chapter  II,  comparison  of  the  experimental 
data  with  those  of  literature  is  presented  in  graphical 
form,  along  with  the  data  obtained  in  this  work  for  regions 
where  no  previous  data  are  available.  Chapter  III,  discusses 
the  theoretical  models  for  prediction  of  equilibrium  partial 
pressures  of  CC^  and  when  other  conditions  are  known. 

This  is  followed  by  comparison  of  the  models. 


CHAPTER  I 


EXPERIMENTAL  METHOD 


In  this  chapter  the  design  considerations,  construc¬ 
tion  and  operation  of  the  experimental  apparatus  used  in 
this  work  will  be  described. 

A.  Design  Considerations: 

Basically  three  different  methods  have  been  employed 
for  vapor  liquid  equilibrium  studies  of  Hydrogen  Sulphide  - 
Carbon  Dioxide  -  Water  -  Ethanolamine  systems: 

(a)  Static  Method 

(b)  Dynamic  Method 

(c)  Circulation  Method 

a)  Static  Method: 

In  this  method  a  closed  cell  containing  the  amine 

solution  with  a  certain  concentration  of  acid  gas  (or  gases) 

is  used.  A  motor  and  gear  train  usually  constitutes  the 

rocking  mechanism.  The  whole  assembly  is  immersed  in  an 

oil  or  air  bath.  This  method  has  been  used  by  Reed  and 
2  8 

Wood  for  measurement  of  CC>2  partial  pressures  ranging 

9 

from  1.5  to  1.7  atm;  Jones  et.al  from  0.56  to  7000  mm  of 

o 

mercury;  Goldman  and  Leibush  from  10  to  2000  mm  of  mercury; 
Lee,  Otto  and  Mather‘S for  partial  pressures  of  C02 
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between  0.1  psia  and  827  psia  and  partial  pressures  of  H2S 
between  0.1  psia  and  306  psia. 

b)  Dynamic  Method: 

Either  a  gas  blend  of  N2 /  CO 2  and  H2S  or  pure  N2  is 
passed  through  a  number  of  saturators  containing  amine  sol¬ 
ution  with  a  certain  concentration  of  acid  gases.  The  gas 
phase  is  repeatedly  analyzed  after  a  certain  time  interval, 
until  two  or  three  similar  readings  indicate  that  equilib¬ 
rium  has  been  reached.  This  method  is  called  the  'Dynamic 
Method' . 

25 

Mason  and  Dodge  used  this  method  to  measure  C02 

partial  pressures  ranging  from  0.19  to  14.9  psia.  Leibush 

23 

and  Shneerson  used  it  for  f^S  partial  pressure  measure¬ 
ments  when  the  partial  pressures  were  less  than  6.7  psia; 
also  for  measuring  solubilities  of  mixtures  of  C02  and  H2S 
when  partial  pressures  of  C02  plus  H2S  were  less  than  3 
psia.  This  method  was  also  used  by  Atwood  et  al?"  when 

partial  pressure  of  H2S  was  between  0.0005  to  5720  mm  of 

24 

mercury.  Lyudkovskaya  and  Leibush  used  it  for  CC>2  par¬ 
tial  pressure  measurements  in  the  range  of  36.8  psia  to 

29 

588  psia:  Riegger,  Tartar  and  Lingafelter  employed  this 
method  for  measurement  of  H2S  partial  pressure  in  the 

26 

range  of  25  to  700  mm  of  mercury.  Muhlbausr  and  Monaghan 

also  used  dynamic  method  for  acid  gas  solubility  studies 
at  2 5°C. 
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c)  Circulation  Method: 

In  this  method  the  vapor  phase  is  circulated  through 
the  equilibrium  cell  in  a  closed  loop.  A  magnetic  pump  or 
some  other  pumping  device  is  used  to  circulate  this  vapor 
phase.  To  analyze  the  vapor  phase  a  small  part  of  the  cir¬ 
culating  vapor  is  bled  to  a  sampling  cell. 

Lee,  Otto  and  Mather  used  this 

method  to  obtain  solubility  data  for  or  CO2  or  their 

mixtures  in  MEA  or  DEA  over  a  wide  range  of  conditions. 

2  6 

Muhlbauer  and  Monaghan  also  used  this  method  for  pressur¬ 
es  ranging  from  1  to  1000  mm  of  mercury. 

27 

Murzin  et  al.  combined  dynamic  and  circulation 

methods  and  termed  it  'Flow-Circulating'  method.  This 

31 

method  is  often  used  in  kinetic  studies  .  In  this  method 
gas  circulates  in  a  closed  system  and  equilibrium  between 
phases  sets  in  during  multiple  bubbling  of  gas  through  the 
fluid  layer.  In  the  meantime,  the  gas  flow  is  continuously 
forced  in  a  closed  cycle  and  equal  flow  equilibrium  gas 
phase  is  withdrawn.  The  gas  flow  forced  in  as  well  as  that 
withdrawn  is  only  a  small  part  (up  to  5%)  of  the  overall 
flow,  therefore  supply  of  pure  gas  into  the  system  does  not 
substantially  alter  the  equilibrium  and  at  the  same  time 
sampling  of  sufficient  quantity  of  gas  is  made  possible. 


. 


' 


7 


B.  Selection  of  a  Method  for  This  Work: 

The  static  method  is  not  suitable  for  measurement  of 

low  partial  pressures  because  the  large  sampling  of  gas 

phase  for  analysis  can  substantially  disturb  equilibrium 

in  the  system.  Besides  at  very  low  partial  pressures  of 

acid  gases  a  very  long  time  is  needed  before  equilibrium 

27 

is  established  .  Since  the  apparatus  in  this  work  had 
to  be  used  for  very  low  partial  pressure  measurement  (of 
acid  gases)  it  was  obvious  at  the  outset  that  this  method 
was  not  reliable. 

The  circulation  method  was  also  rejected  because  for 
measurement  of  small  partial  pressures  it  requires  the  cir¬ 
culation  of  a  large  quantity  of  gas  and  this  makes  the 

2  6  27 

method  excessively  cumbersome  ' 

The  dynamic  method  has  the  advantage  that  a  suffic¬ 
ient  quantity  of  the  vapor  phase  is  continuously  available 
for  analysis.  Thus  gas  phase  analysis  in  no  way  affects 
the  equilibrium.  Also,  if  the  flow  rate  of  the  carrier  gas 
(e.g.  N2  in  this  case)  is  chosen  carefully  equilibrium  is 
reached  within  a  reasonable  time  interval.  This  method 
has  the  disadvantage  that  at  higher  temperatures ,  the  amine 
concentration  changes  due  to  evaporation  of  water  unless  a 
high  gauge  pressure  is  maintained.  But  the  amine  concent¬ 
ration  changes  can  be  controlled  to  within  reasonable  limits 


,  .* 


by  a  slightly  higher  absolute  pressure  than  the  partial 
pressure  of  water  plus  amine  at  that  temperature  and  by 
maintaining  a  low  flow  rate  of  ^  through  the  saturators. 

Two  alternate  routes  are  available  in  the  dynamic 
method  itself: 

a)  Solution  Saturation  Method: 

Here  a  carrier  gas  of  known  composition  viz.  ^ 
containing  acid  gas(es)  is  passed  through  a  series  of  sat¬ 
urators.  Equilibrium  is  indicated  when  inlet  and  outlet 
gas  compositions  are  identical. 

b)  Gas  Saturation  Method: 

Pure  carrier  gas  is  passed  through  a  number  of 

saturators  containing  amine  solution  with  acid  gas(es). 

Since  this  work  had  to  be  carried  out  in  the  region 
where  partial  pressures  of  CC^  and  H2S  were  always  less 
than  1  psia  the  gas  saturation  method  was  chosen.  The  ob¬ 
vious  reason  being  that  it  is  not  practical  to  have  a 
uniform  concentration  of  acid  gas(es)  in  the  carrier  gas 
when  the  partial  pressure  of  these  acid  gases  is  extremely 
low^  (e.g.  below  0.1  psia). 
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C .  Description  of  the  Apparatus: 

The  apparatus  used  for  this  work  was  quite  simple 
and  is  shown  schematically  in  Figure  1.  Basically  it 
consisted  of  three  identical  250  ml  gas  washing  bottles 
connected  in  series.  These  bottles  were  kept  in  an  oil 
bath.  Carrier  gas  (N2)  passed  through  a  25  foot  long, 

0.25  in.  O.D.  stainless  steel  coil  immersed  in  the  oil 

along  with  most  connections  and  all  the  three  bottles. 

After  passing  through  the  heating  coil  N2  was  bubbled 
through  the  three  saturators  in  series.  N2  flow  was  regul¬ 
ated  by  a  needle  valve  to  make  sure  that  the  flow  rate  was 
constant.  The  three  saturators  contained  aqueous  amine 
solution  with  a  certain  concentration  of  acid  gas(es). 

This  solution  was  prepared  by  passing  the  acid  gas(es) 
through  the  amine  solution  in  a  separate  bottle.  After  the 
saturators  N2  saturated  with  water  and  acid  gas(es)  passed 
through  a  trap  to  remove  water  droplets  before  going  to  the 
gas  chromatograph  (GC) ,  and  finally  to  the  exhaust  line. 

The  gas  continuously  passed  through  the  GC  sampling  loops; 
therefore  gas  phase  monitoring  was  possible  at  whatever  time 
intervals  desired,  without  disturbing  the  system.  The  pres¬ 
sure  in  the  final  saturator  was  indicated  by  a  mercury  or 
water  manometer  depending  on  whether  high  or  low  gauge  pres¬ 
sure  was  being  maintained.  Liquid  samples  were  withdrawn 
from  the  last  bottle  by  a  back- forcing  method  to  be  describ¬ 


ed  later. 
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All  interconnecting  lines  were  of  seamless  stainless 
steel  tubing  except  at  inlet  and  outlet  of  the  gas  washing 
bottles  which  had  glass  ball  and  socket  joints.  Glass  and 
metal  were  connected  by  Kovar-to-borosilicate  glass 
seals.  Only  stainless  steel  valves  and  connections 
were  used  because  of  the  corrosive  nature  of  the  acid  gases 
(especially  I^S)  and  the  amine  solutions. 

Glass  tubing  fused  from  the  side  of  the  bottle  and 
reaching  the  bottom  of  the  bottle  (from  inside)  was 
used  for  filling  and  discharging  of  the  amine  solutions. 
These  lines,  too,  were  connected  by  Kovar-to-borosilicate 
glass  seals  to  stainless  steel  tubing  through  a  ball  and 
socket  joint.  The  operation  of  filling,  discharging  and 
cleaning  of  the  bottles  without  removing  the  apparatus  from 
the  oil-bath  will  be  described  later  in  the  section  on  ex¬ 
perimental  procedure. 

The  traps  used  before  the  GC  and  the  manometer  to 
remove  condensed  water  drops  were  simply  250  ml  erlenmeyer 
flasks  with  a  side  arm  which  served  as  an  outlet  line.  The 
inlet  stainless  steel  tubing  passed  through  an  air  tight 
rubber  stopper.  The  outlet  lines  from  the  traps  were  of 
polyethylene  plastic  because  the  vapor  phase  had  cooled 
to  near  room  temperature.  The  trap  near  the  manometer  was 
used  to  prevent  any  vapor  from  reaching  the  column  of  mer¬ 
cury  if  a  leak  developed  in  the  line  connecting  the  last 
bottle  to  the  manometer. 


. 


For  indication  and  approximate  measurement  of  gas 

flow  rate  through  the  system,  a  rotameter  was  used  before 
the  carrier  gas  entered  the  saturators.  A  similar  rota¬ 
meter  was  initially  used  before  the  gas  entered  the  GC, 
but  this  had  to  be  replaced  by  a  soap-bubble  flow  meter, 
because  of  the  moist  nature  of  the  gas  which  caused  the 
indicator  to  stick  in  place  at  low  flow  rates.  A  flow 
meter  is  desirable  after  the  GC  because  in  the  event  of  a 
leak  in  the  apparatus  the  two  flow  meters  (one  before  the 
saturators  and  the  other  after  them)  would  give  different 
flow  rates.  The  soap-bubble  flow  meter  (after  the  GC)  was 
found  to  be  adequate.  The  only  disadvantage  in  using  a 
soap-bubble  flow  meter  was  that  it  sometimes  affected  the 
pressure  in  the  system  since  the  line  to  the  exhaust  had 
to  be  closed  to  divert  flow  through  the  bubble-meter . 

The  saturators  were  mounted  on  a  stainless  steel 
platform  with  a  heating  coil  surrounding  them  on  supports 
from  the  platform.  This  assembly  was  immersed  in  oil  but 
all  the  valves  were  above  the  oil  surface  and  could  be 
reached  easily.  Filling-discharging  lines  and  the  liquid 

sampling  lines  could  also  be  reached  conveniently. 

% 

The  bath  consisted  of  a  3  foot  x  2  foot  x  1.5  foot 
rectangular  stainless  steel  container  containing  about  40 
gallons  of  Grade-30  motor  (mineral)  oil  as  the  heating 
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medium.  The  outer  walls  were  covered  by  styrofoam  insula¬ 
tion  and  the  bath  was  enclosed  in  a  hardboard  casing.  Agi¬ 
tation  of  oil  was  obtained  by  a  stirrer  rotated  by  a  0.25 
HP  electric  motor  mounted  on  the  cover  of  the  bath.  Heat¬ 
ing  was  accomplished  by  a  1500  watt,  220  volt  cylindrical 
heating  element  (General  Electric)  which  was  screwed 
into  the  bath  about  4  inches  above  the  bath  floor.  The 
heater  was  connected  to  a  Hallikainen  Instrument  Thermo- 
trol,  which  controlled  the  bath  temperature  to  within 
±0.01°C.  The  sensor  from  the  temperature  controller  was 
passed  into  the  oil  through  a  hole  in  the  bath  cover.  The 
temperature  of  the  bath  was  indicated  by  two  pre-calibrated 
thermometers , one  in  the  middle  and  the  other  at  the  end  of 
the  cover. 

The  gas  inlet,  gas  outlet  and  the  line  to  the  mercury 
(or  water)  manometer  came  out  of  the  bath  through  two  holes 
in  the  bath  cover.  These  lines  had  stainless  steel  quiek- 
conneots  so  that  the  apparatus  inside  the  bath  could  be 
disconnected  conveniently  whenever  needed.  A  small  lid 
just  above  the  saturator-assembly  allowed  access  to  the 
valves  and  for  liquid  sampling,  filling,  discharging  and 
cleaning  of  saturators.  At  other  times  the  lid  was  kept 
closed  to  minimize  the  escape  of  oil  vapors. 


»f  6 
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D.  Operating  Details: 

a)  Start-Up: 

The  bottles  were  filled  with  desired  amine  solution 
and  lowered  on  the  stainless  steel  platform  which  held  them 
in  place.  The  ball  and  socket  joints  were  properly  greased 
(with  Dow  Corning  high  vacuum  silicone  grease)  and 
clamped.  The  whole  assembly  was  lowered  into  the  oil  bath 
with  only  the  inlet,  outlet  and  manometer  lines  coming  out 
of  the  top  of  the  bath  through  two  holes.  These  lines  were 
connected  to  the  manometer,  inlet  gas  line  and  outlet  line 
to  the  GC,  through  quick-connects.  At  this  time  all  the 
valves  were  in  the  closed  position.  The  stirrer  was  turned 
on  and  after  some  turbulence  in  the  oil  had  been  obtained, 
the  heater  along  with  the  thermotrol  (temperature  control¬ 
ler)  were  turned  on.  The  thermotrol  reading  was  set 
through  a  calibration  curve  obtained  earlier.  It  took  a 
few  hours  for  the  bath  to  reach  a  desired  temperature.  For 
example,  it  took  about  6  to  7  hours  for  the  bath  to  reach 
100 °C  from  room  temperature. 

b)  Establishment  of  equilibrium: 

Once  the  temperature  of  the  bath  had  reached  a  steady 

value,  carrier  gas  (^J  flow  through  the  saturators  was 
started.  A  nitrogen  pressure  of  about  5  psig  was  first 
set.  Then  valves  1,  2,  3,  4,  5,  6  and  7  were  opened  in 


■ 


15 


that  order.  A  small  change  in  mercury  level  in  the  manometer, 
M-2,  showed  that  gas  had  started  flowing  through  the  system. 
Also  the  rotameter  showed  flow  of  N2  through  the  system. 

Now  valve-13  was  opened  and  then  valve-7  closed,  to  measure 
the  flow  rate  of  carrier  gas  through  the  saturators  by  the 
soap-bubble  meter.  If  the  flow  rate  was  above  15  ml/min 
valve-1  was  closed  slightly  until  the  rate  reached  about 
this  value.  The  final  flow  rate  setting  was  obtained  by  the 
fine-metering  valve- 6.  A  N2  flow  rate  of  about  10  ml/min 
was  found  to  be  the  optimum.  Once  the  desired  flow  rate 
had  been  obtained  valve-7  was  opened  again  and  valve-13 
closed.  Valve-8  was  now  opened  to  connect  the  last  satura¬ 
tor,  B-3,  to  the  manometer,  M-l,  which  gave  the  gauge  pres¬ 
sure  in  B-3. 


For  temperatures  up  to  80°  C,  N2  flow  through  the 

amine  solution  did  not  affect  the  amine  concentration  much 

because  the  partial  pressure  of  2.5  normal  aqueous  MEA  sol- 

7  7 

ution  is  about  2.8  psia  at  60°C  and  6.6  psia  at  80°C. 
About  4  hours  were  allowed  for  the  N2  to  flow  through 

the  saturators.  After  that  the  gas  composition  was 
monitored  every  half-hour  using  the  GC.  The  gauge  pressure 
from  M-l  was  also  noted  each  time.  When  two  identical  GC 
analysis  had  been  obtained  and  M-l  readings  were  also  al¬ 
most  the  same  equilibrium  had  been  established.  To  make 
sure  one  hour  more  was  allowed  and  M-l  and  GC  readings 
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taken  again.  If  these  readings  were  similar  to  the  previous 
ones  establishment  of  equilibrium  was  confirmed.  Liquid 
samples  were  now  taken  using  the  method  described  later. 

For  a  temperature  of  100°C  the  partial  pressure  of 

7 

2.5  normal  aqueous  MEA  solution  was  about  13.7  psia  which 
was  slightly  above  atmospheric  pressure.  Therefore  to 
avoid  excessive  removal  of  water  and  hence  appreciable 
change  in  amine  concentration  an  absolute  pressure  of  about 
15  to  16  psia  was  maintained  in  B-3.  This  was  achieved  by 
closing  valve-5  slightly  in  combination  with  valve-6  to 
maintain  a  flow  rate  of  carrier  gas  around  10  ml/min.  In 
the  case  of  5  normal  MEA  solution  the  partial  pressure  of 

7 

water  plus  MEA  was  about  13  psia  .  Therefore,  an  absolute 
pressure  of  about  14.5  to  15  psia  kept  the  amine  normality 
change  within  desirable  limits.  N«  was  allowed  to  flow 
through  the  saturators  at  10  ml/min  for  about  3  hours  be¬ 
fore  starting  the  monitoring  of  the  vapor  phase.  Monitor¬ 
ing  was  continued  every  half  hour  until  consistent  readings 
were  obtained.  A  close  watch  was  kept  on  M-l  to  make  sure 
that  gauge  pressure  did  not  vary  much  between  monitorings 
because  a  change  in  pressure  disturbs  the  equilibrium. 
Micrometering  valve-6  was  used  to  maintain  this  pressure 
to  a  near  constant  value.  Once  two  consistent  readings  of 
GC  and  M-l  were  obtained,  equilibrium  had  been  reached  and 
liquid  sample  was  now  taken.  The  time  needed  to  reach 
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equilibrium  depends  on  the  acid  gas  concentration  in  the 
liquid  phase.  For  higher  concentrations  it  usually  took 
about  4  to  5  hours  if  the  gauge  pressure  in  B-3  was  main¬ 
tained  around  a  certain  value.  For  lower  acid  gas  concen¬ 
trations  a  little  longer  time  was  usually  needed. 

c)  Liquid  sampling: 

Once  the  GC  analysis  of  the  vapor  phase  had  shown 
that  equilibrium  had  been  reached  and  also  a  reading  of  the 
manometer,  M-l,  had  been  taken,  it  was  time  for  obtaining 
a  representative  sample  of  the  liquid  phase.  Saturator  B-3 
was  the  one  with  which  the  gas  phase  was  in  contact  last 
and  hence  the  liquid  sample  was  withdrawn  from  B-3.  To  do 
this  valves  2,  4,  5  and  8  were  closed.  Valve-9  was  now 
opened  to  apply  a  slight  back  pressure  on  the  liquid  in  B-3. 
Valve-10  was  now  opened.  Liquid  was  forced  through  the 
fritted  glass  end  in  B-3  through  valve-10  and  into  a  liquid 
sampling  bottle  at  port,  p-10.  The  first  5  to  10  ml  of  the 
sample  were  discarded  since  they  had  flushed  the  line  and 
contained  some  drops  from  the  previous  sample.  After  this 
flushing  about  10  ml  were  collected  in  a  liquid  sampling 
bottle  which  was  then  stoppered  and  allowed  to  cool  to  room 
temperature  before  chemical  analysis. 

d)  Discharging,  Cleaning  and  Refilling  of  Saturators: 

To  discharge  amine  solution  the  bottles  were  emptied 
in  the  order:  B-3,  B-2  and  finally  B-l.  Valve-9  and  then 


. 
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the  three-way  valve-12  for  B-3  were  opened.  All  the  other 
valves  except  1  were  closed  earlier.  The  discharged  solu¬ 
tion  was  collected  in  the  filling-bottle  or  a  container  at 
port,  p-12.  Once  B-3  had  been  emptied  valve-4  was  opened 
and  valve-12  was  turned  towards  B-2.  Contents  of  B-2  were 
again  collected  at  port  p-12.  After  B-2  had  been  emptied 
valves  3  and  then  11  were  opened  to  discharge  the  contents  of 
B-l  at  port,  p-11. 

Filling  of  the  bottles  could  be  accomplished  in  any 
order  desired.  To  fill  B-l,  a  filling-bottle  was  connected 

f 

to  p-11.  Valves  11,  3,  4,  5,  6  and  7  were  then  opened  in 
this  order.  Amine  solution  (or  cleaning  liquid)  was  forced 
into  B-l  by  applying  Nj  pressure  above  the  liquid  in  the 
filling-bottle .  Once  B-l  had  been  filled  with  a  desired 
solution,  valves  11  and  3  were  closed.  The  filling-bottle 
was  now  connected  to  p-12  and  valve-12  was  opened  for  B-2. 

N2  pressure  on  liquid,  in  the  filling-bottle ,  forced  it 
into  B-2.  After  B-2  had  been  charged  valve-4  was  closed 
and  valve-12  opened  for  B-3.  B-3  was  filled  in  the  same 

way  as  the  other  two  bottles. 

For  cleaning  of  the  bottles  warm  distilled  water 
(or  amine  solution)  was  pushed  into  each  bottle  using  the 
technique  described  for  filling  of  the  bottles.  N2  was 
bubbled  through  this  liquid  for  a  few  minutes  (for  proper 
rinsing)  and  then  the  liquid  was  removed  using  the  proce¬ 
dure  described  for  discharging. 


■ 

. 


19 


Note  that  when  the  saturator  assembly  was  at  a  high 
temperature  (above  80°C)  cleaning  liquid  or  new  amine  sol¬ 
ution  was  warmed  before  charging  to  prevent  cracking  (or 
breaking)  of  the  glass  due  to  sudden  contraction. 

E.  Analytical  Methods : 

a)  Gas  Phase  Analysis: 

A  model  5710A  Hewlett-Packard  gas  chromatograph 
equipped  with  a  Model  1750A  disc  integrator  recorder  was 
used  for  all  gas  phase  analyses.  The  column  was  a  10 
foot  long,  0.25  in.  O.D.  stainless  steel  tube  filled  with 
CHROMOSORB-104  packing.  The  column  was  kept  at  130°C 
for  H2S  and  110 °C  for  CC>2  or  for  CC>2  and  H2S  mixtures. 

b)  Liquid  Phase  Analysis: 

In  this  work  both  direct  titration  methods  and 
methods  which  rely  on  the  evolution  of  C02  and/or  H2S 
were  tested  for  the  chemical  analysis  of  the  acid  gases. 

It  was  found  that  for  very  small  concentrations  (less 
than  0.1  mole  acid  gas(es)  per  mole  of  amine)  of  C02  and 
H2S  the  evolution  methods  are  not  reliable  since  the 
quantity  of  sample  available  for  analysis  was  small 
(about  10  ml) .  Also  almost  always  a  fraction  of  the  acid 
gas(es)  remained  in  the  liquid  no  matter  how  much  the 
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technique  was  refined.  A  brief  description  of  the  analysis 
methods  tried  follows :- 

Method-1 : 

Dow  Method;  Iodometric  Determination  of  evolved  H2S: 
H2S  was  evolved  from  an  aliquot  portion  of  the  liquid 
sample  by  addition  of  concentrated  sulphuric  acid  and  boil¬ 
ing.  Purified  air  was  drawn  in  by  aspiration  through  the 
evolution  flask  carrying  the  evolved  H2S  through  a  reflux 
condenser,  then  through  a  gas  scrubber  containing  a  known 
volume  of  standard  iodine  solution.  Excess  iodine  was 
determined  by  titration  with  standard  sodium  thiosulphate 
solution. 

Method- 2 : 

H2S  by  Iodine-Thiosulphate  Titration:-  The  method 

14 

is  that  outlined  by  Kolthoff  and  Sandell  .  An  aliquot 
portion  of  the  liquid  sample  was  put  into  acidified  iodine 
solution.  Excess  iodine  was  back  titrated  using  standard 
sodium  thiosulphate  solution.  Sulphuric  acid  in  iodine 
should  be  at  least  18  moles  per  mole  of  DEA  or  11  moles 
per  mole  of  MEA. 

Method- 3 : 

2 

H2S  by  Arsenite- Iodine  Titration  :-  Here  Arsenious 
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sulphide  was  precipitated  by  H2S  and  the  excess  arsenite 
was  back  titrated  using  iodine  solution.  The  reaction  for 
this  precipitation  is: 

As203  +  3H2S  — ♦  as2S3  +  3H2° 

The  problem  with  this  method  was  that  a  large  volume 
(about  20  ml)  of  sample  was  needed  for  reliable  results. 
With  about  2  to  3  ml  liquid  samples  reproducible  results 
were  not  obtained. 

Method- 4 : 

H2S  by  direct  I2  titration:  An  aliquot  of  sample 
was  titrated  directly  with  0.1  normal  standard  iodine 
solution  in  a  medium  buffered  with  sodium  bicarbonate. 

Method- 5 : 

g 

Dow  method  for  CC>2  :  C02  was  evolved  from  an  ali¬ 

quot  portion  of  the  sample  by  addition  of  concentrated 
sulphuric  acid  and  heating.  Evolved  C02  was  carried 
through  a  reflux  condenser  by  purified  air  drawn  in 
through  an  aspirator  and  then  absorbed  in  a  known  volume 
of  0.1  N  Barium  Hydroxide  (Ba(OH)2).  The  excess  Ba(OH)2 
was  determined  by  titration  with  0.1  normal  hydrochloric 
acid  to  the  phenolphthalein  end  point.  Again  the  problem 
with  this  method  was  that  to  use  up  a  reasonable  amount  of 
Ba(OH)2  a  large  amount  of  liquid  sample  was  needed. 


> 

' 
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especially  for  very  low  concentrations  of  CC>2.  Since  a 
large  liquid  sample  was  not  feasible,  this  method  was  dis¬ 
carded. 

Method- 6 : 

CC>2  by  precipitation  as  BaCC>3 :  The  CO 2  in  an 
aliquot  portion  of  the  sample  was  precipitated  as  Barium 
Carbonate  (BaCO^)  in  the  presence  of  excess  Barium  Chloride 
(BaCl2) .  The  solution  was  filtered  and  the  BaCO^  precipi¬ 
tate  was  titrated  with  0.1  normal  hydrochloric  acid  using 
modified  methyl  orange  (Methyl  Orange  -  Xylene  Cyanol 
solution)  indicator. 

Determination  of  Amine  Concentration:  The  amine 
concentration  in  the  liquid  was  determined  by  titration  of 
an  aliquot  of  the  sample  with  about  0.1  normal  sulphuric 
acid  solution  employing  methyl  red  as  an  indicator. 

F.  Selection  of  Methods  for  Liquid  Analysis: 

As  mentioned  earlier  methods  relying  on  evolution 
of  acid  gas(es)  were  not  reliable  at  very  low  concentra¬ 
tions  of  acid  gases  in  the  amine  solution.  Therefore, 
these  methods  were  not  pursued  after  preliminary  tests. 

Out  of  the  direct  titration  methods  for  H2S,  Method-4  al¬ 
ways  gave  a  larger  value  of  H2S  composition  in  the  amine 
solution  as  compared  to  iodine-sodium- thiosulphate  method 
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(Method- 2) .  Data  in  the  literature  were  reproduced  within 
reasonable  limits  using  the  iodine- thiosulphate  method, while 
with  Method-4  the  f^S  partial  pressure  was  always  lower 
for  the  determined  H2S  composition  using  this  method.  This 
indicated  that  liquid  phase  determination  was  incor¬ 

rect  in  the  positive  direction  i.e.,  Method-4  always  gave  a 
higher  H2S  composition  than  actual.  Later  it  was  found 

that  this  was  due  to  the  fact  that  some  iodine  was  used  up 

2- 

by  thiosulphate  (S2C>3  )  present  in  a  sample.  If  the  amount 

of  iodine  used  up  by  thiosulphate  was  determined  and  then 

subtracted  from  the  total,  the  iodine  consumed  by  H2S  was 

obtained.  This  gave  results  which  were  nearly  identical 

to  those  given  by  iodine- thiosulphate  (Method- 2)  method. 

Since  reliability  of  iodine-thiosulphate  titration  has  been 

14 

confirmed  before  and  since  determination  of  thiosulphate 
in  a  sample  is  quite  cumbersome,  the  iodine-thiosulphate 
method  (Method- 2)  was  adopted  and  used  for  all  the  H2S 
analysis . 

As  for  CC>2  Method- 6  was  found  to  be  good.  It 
reproduced  data  in  the  literature  reasonably  well  and  also 
gave  reproducible  results.  Its  only  disadvantage  is  that 
it  is  slightly  cumbersome.  Because  of  absence  of  any 
other  reliable  and  less  cumbersome  technique  for  CC>2 
analysis  in  the  literature,  this  method  was  used  for  all 
C02  determinations. 


The  method  used  for  amine  analysis  is  the  one  used 
by  Girdler  Corporation,  Louisville,  Kentucky  (SM-1M  for 
MEA;  SM-lD  for  DEA) .  This  method  was  tested  by  preparing 
standard  amine  solution  and  was  found  to  be  very  accurate 
and  reliable. 


G.  Details  of  the  Methods  Used  for  CC>2,  H2S 
And  Amine  Determinations: 

a)  CC>2  Determination  (Method- 6) 

Procedure:  About  20  to  25  ml  of  distilled  water 

was  heated  in  a  250  ml  erlenmeyer  flask.  Once  the  water 
had  started  boiling  approximately  2  grams  of  solid  BaCl2 
were  added.  The  flask  was  swirled  to  dissolve  the  BaCl2. 
Immediately  by  a  syringe  2  ml  of  the  sample  solution  was 
added  to  the  flask.  The  flask  was  removed  from  heat  and 
allowed  to  stand  for  about  five  minutes  for  the  white  pre¬ 
cipitate  to  settle  down  somewhat.  The  flask  was  then 
stoppered  and  cooled  under  cold  tap  water  for  about  ten 
minutes.  Three  7  cm  41  Whatman  ashless  filter  papers 
were  placed  in  series  in  a  7  cm  Buchner  funnel  placed  in  a 
large  erlenmeyer  flask  with  the  sidearm  attached  to  a  water 
aspirator.  Slowly  the  precipitate  was  poured  on  to  the 
filter  papers.  The  flask  containing  the  precipitate  was 
washed  with  two  25  ml  portions  of  distilled  water.  Each 
washing  was  transferred  to  the  filter  paper.  The  filtrate 
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was  discarded.  The  filter  papers  with  the  washed  precipi¬ 
tate  were  transferred  to  the  erlenmeyer  flask  which  had 
the  precipitate  previously  and  about  fifty  ml  of  distilled 
water  were  added  to  it.  The  contents  of  the  flask  were 
swirled  with  a  magnetic  stirrer,  until  the  filter  paper 
had  been  torn  into  small  pieces  (this  puts  almost  all  the 
precipitate  in  solution  and  therefore,  it  is  easier  to  tit¬ 
rate)  .  Approximately  five  drops  of  modified  methyl  orange 
indicator  (Methyl  Orange  -  Xylene  Cyanol  solution)  were 
added  to  the  flask  and  the  mixture  titrated  with  about  0.1 
normal  hydrochloric  acid.  The  mixture  was  titrated  to  a 
grey  end  point.  Any  barium  carbonate  precipitate  sticking 
to  the  sides  of  the  flask  was  washed  down  with  distilled 
water  and  the  titration  continued  to  the  same  grey  end 
point. 

b)  I^S  Determination  (Method-2) 

Procedure:  20  ml  of  0.1  normal  standard  iodine 

solution  was  pipetted  into  a  500  ml  erlenmeyer  flask.  The 
proper  excess  of  sulphuric  acid  for  the  size  of  amine 
sample  aliquot  was  added.  (The  ratio  of  sulphuric  acid  to 
DEA  should  be  at  least  18  moles  I^SO^  per  mole  of  DEA  and 
11  moles  sulphuric  acid  per  mole  of  MEA) .  The  magnetic 
stirring  bar  was  now  put  in  and  by  a  syringe  2  ml  aliquot 
of  sample  was  added  while  stirring.  The  excess  iodine  was 
back  titrated  with  standard  0.1  normal  sodium  thiosulphate 
solution  to  a  light  yellow  colour.  About  5  ml  of  0.2 


. 
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percent  starch  indicator  solution  was  now  added  and  the 
titration  continued  to  the  change  from  blue-black  to 
water-white  colour. 

c)  MEA  (or  DEA)  Determination: 

Procedure:  To  an  erlenmeyer  flask  containing  about 

100  ml  of  distilled  water,  two  ml  of  sample  solution  were 
added.  A  few  drops  of  methyl  red  indicator  were  added  un¬ 
til  the  colour  was  distinctly  yellow.  The  solution  was 
then  titrated  using  ca.  1  normal  sulphuric  acid.  The  end 
point  was  indicated  by  change  in  colour  of  solution  from 
yellow  to  pink. 


CHAPTER  II 


EXPERIMENTAL  RESULTS  AND  THEIR  COMPARISON 

This  chapter  presents  in  graphical  form,  the  exper- 
mental  data  obtained  in  this  work. 

To  test  the  apparatus,  initial  measurements  were 

made  for  conditions  of  solution  strength  and  temperatures 

for  which  some  data  were  available  in  the  literature. 

Measurements  were  made  for  pure  H2S  and  pure  C02  in  2.5  N 

MEA  at  80°  and  100°C.  For  the  two  100°C  isotherms  and  the 

80°C  isotherm  for  H2S,  some  data  points,  for  low  loadings 

of  acid  gases  (ca  0.1  moles/mole  of  MEA  and  less),  have 

9  33 

been  obtained  by  Jones  et  al.  ,  and  Garst  and  Lawson 
These  points  have  been  plotted  on  the  same  figures.  In  this 
work,  data  were  also  obtained  for  mixtures  of  CC>2  and  H2S 
in  5  N  MEA  at  100°C.  The  loadings  of  H2S  and  CC>2  were  less 
than  ca.  0.12  moles/mole  of  MEA.  For  these  low  acid  gas 
concentrations,  no  data  have  been  published  to  date. 

Figures  2  and  3  present  the  80 °C  isotherms  for 
pure  H2S  and  CC>2  in  2.5  N  MEA.  Figures  4  and  5  show  the 
100 °C  isotherms  for  pure  H2S  and  CC>2  in  2.5  N  MEA.  Smooth¬ 
ed  mixture  data,  for  CC>2  and  H2S  in  5  N  MEA  at  100°C,  are 
plotted  in  Figures  6  and  7.  The  raw  experimental  data  are 
tabulated  in  Appendix  B. 
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Figure  2  -  Solubility  of  H2S  in  2.5  N  MEA  Solution 
at  80 °C . 
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Solubility  of  C02  in  2.5  N  MEA 
Solution  at  80°C. 
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Solubility  of  I^S  in  2.5  N  MEA 
Solution  at  100°C. 


Figure  4 
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Figure  5  -  Solubility  of  C02  in  2.5  N  MEA  Solution 


at  100°C 
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Figure  6  - 


H2S 

Effect  of  C02  on  the  Partial  Pressure 
of  H2S  Over  5  N  MEA  Solution  at  100 °C, 
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Effect  of  on  the  Partial  Pressure 
of  CC>2  Over  5  N  MEA  Solution  at  100°C. 


Figure  7 


CHAPTER  III 


THEORETICAL  PREDICTION  OF  PARTIAL  PRESSURES  OF  H2S  and  C02 

While  the  solubility  data  are  useful  in  themselves, 
it  is  desirable  to  have  a  computer-oriented  correlation 
which  reproduces  the  experimental  results.  It  would  also  be 
useful  for  regions  where  no  data  are  presently  available. 

A  number  of  researchers  have  attempted  to  correlate 
the  solubility  data  for  C02  and/or  H2S  in  ethanolamine  sol¬ 
utions.  A  common  starting  point  for  all  the  thermodynamic 
models  proposed  so  far  is  the  postulation  of  reactions  oc- 
curing  in  the  solution.  For  simple  systems,  i.e.  systems 
containing  either  H2S  or  CC>2  in  aqueous  ethanolamine  (EA) 
solution,  correlations  have  been  proposed  by  different 
authors.  Atwood  et  al.^  proposed  a  correlation  for  pre¬ 
dicting  partial  pressures  of  H2S  in  aqueous  ethanolamine 
solutions  based  on  chemical  reactions  taking  place  in  the 
solution.  Klyamer^  made  slight  modifications  in  the 

above  correlation  and  has  also  made  some  comparisons  with 

12 

experimental  data.  Klyamer  and  Kolesnikova  discussed  a 
general  mathematical  description  of  experimental  data  for 
the  thermodynamic  equilibrium  in  the  C02-water-MEA (DEA) 
systems.  In  postulating  the  equilibrium  equations  they 
have  accepted  the  assumptions  made  by  Atwood  et  al.^. 

They  found  that  the  maximum  error  in  the  calculated  values 
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. 


; 


compared  to  the  experimental  ones  for  a  given  set  of  condi 
tions  was  not  larger  than  20  percent  of  CC>2  partial  pres¬ 
sure.  Theoretical  treatment  of  equilibria  for  the  absorp¬ 
tion  of  CC>2  into  aqueous  ethanolamine  solutions  has  also 

4 

been  discussed  by  Danckwerts  and  McNeil  . 

13 

More  recently  Klyamer  et  al.  have  proposed  a 
thermodynamic  model  for  predicting  equilibrium  partial 
pressures  of  acid  gases  in  aqueous  ethanolamine  solutions. 
They  have  obtained  equations  which  describe  the  functional 
dependence  between  the  partial  pressure  of  the  acid  gas 
components  and  the  composition  of  the  solution  at  various 
temperatures.  In  essence,  their  model  combines  the  chem¬ 
ical  equations  postulated  for  the  simple  systems,  I^S- 
water-EA  and  CC^-water-EA,  and  also  incorporates  the 
assumptions  made  for  these  simple  systems  by  other  authors 

Most  recently  (while  this  work  was  in  progress) 

Kent  and  Eisenberg^^  presented  a  modified  version  of  the 

13 

Klyamer  et  al.  model.  This  model  is  essentially  based 
on  postulation  of  similar  reactions  as  those  of  the 
Klyamer  et  al.  model;  however,  all  the  non-idealities 
have  been  incorporated  into  only  two  parameters. 

13 

Both  the  Klyamer  et  al.  and  the  Kent  and 
Eisenberg"^  Models  are  discussed  in  detail  in  the  follow¬ 
ing  pages  along  with  the  modification  in  Klyamer  et  al. 


model . 


. 
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A.  Correlation  Based  on  the  Thermodynamic  Model  of 

13 

Klyamer  et  al.  ; 

The  assumptions  made  in  postulating  the  chemical 
reactions  are  listed  below: 

a)  The  salts  (bisulphides,  carbonates  and  carbamates) 
which  are  created  by  the  absorption  of  H2S  and  CC>2  are 
fully  dissociated'*' '  ^  ^  # 


b)  The  activity  coefficients  of  different  ions  in  sol¬ 
ution  are  equal  and  independent  of  temperature'*'. 

c)  The  activity  of  water  is  equal  to  its  molarity'*'. 


d)  The  physical  solubility  of  C02  or  H2S  is  not 
changed  by  the  addition  of  ethanolamine  (EA)  and  salts  of 


EA 


1,12 


On  the  basis  of  the  above  assumptions  the  follow¬ 
ing  reactions  for  the  I^S-CO^water-MEA  systems  can  be 
written: 


RNH2  +  H20  RNH3+  +  OH 
H2S  ^=±  H+  +  HS” 

h2o  ^=±  n+  +  oh” 


(a) 

(b) 

(c) 


HS 


(d) 


. 

. 
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2RNH2  +  C02  ^  RNH*  +  RNHCOO  - (e) 

C02  +  H20  ^=±  E+  +  HCO~  - (f) 

HCO”  H+  +  C032”  - (g) 

where  R  stands  for  -CH2-CH2~OH  group.  Analogous  equations 
may  be  written  (to  represent  the  reactions)  for  the 
H2S-C02-water-DEA  system. 

The  following  balance  equations  for  the  reacting  compounds 
can  be  formed: 


Charge  balance  (or  Electroneutrality)  equation: 

[RNH3]+[H+]  =  [HC03]+ [RNHC00]+2 [C032"]+ [0H”]+ [HS~]+2 [S2"] 

- (!) 


Mass  balance  equations: 


[RNH2]  +  [RNHCOO]  +  [RNHj]  =  m  - (2) 

[HS~]  +  [H2S]  +  [S2- ]  =  mXH  s  - (3) 

2 

[HCO~]  +  [C02]  +  [C032-]  +  [RNHCOO]  =  mXCQ  - (4) 

where. 


m  =  moles  of  EA/1000  gm  of  water 

Xc0  =  moles  of  C02/mole  of  EA 
2 

X„  c  =  moles  of  H0S/mole  of  EA 
H  ^  o  Z 


' 


' 
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[Ci]  =  moles  of  compound  (or  ion)  C^/1000  gm  of 
water  (Molality) . 


The  following  thermodynamic  equilibrium  constant  expressions 
can  also  be  written: 


2  [RNH*]  [0H_] 
_  Y _  ^ 

•  r\  ~ 

1  aa  [RNH21 


=  Y2  [H+l  [OH  ] 

a 

=  Y2  [H~*~  ]  [HS~] 

[h2s] 

=  Y  tH+]  [S2~] 

[HS_] 

=  [H2S1 
PH25 

2  [RNH*]  [RNHCOO] 
=  JL  •  _ 

a2a2  [RNH2]2  PC02 

Y2  [H+]  [HCO“] 
a  [C02] 

y  [H+]  [C032~] 

[HCO~] 


[co2] 


-(5) 

-(6) 

-(7) 

-(8) 

-(9) 

(10) 


(ID 

(12) 


(13) 
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WHERE 

[RNHCOO],  [HCOj],  [C032_],  [C02],  [H2S],  [RNH+] , 

[RNH2],  [OH  ],  [H+] ,  [HS  ]  and  [S2  ]  are  respectively,  the 

concentrations  (molalities)  of  the  carbamate,  bicarbonate, 
carbonate  ions,  the  physically  dissolved  CO 2  and  H2S,  con¬ 
centrations  of  ionized  and  un-ionized  ethanolamine ,  hydroxyl 
ions,  hydrogen  ions,  bisulphide  and  sulphide  ions.  PCQ 

2 

and  e  are  the  partial  pressures  of  C0~  and  H0S  respect- 

11 2  ^  m  Z 

ively  in  mm  of  mercury;  Km,  Kly,  K2y,  Kir  K^,  H^ ,  HH2S' 

K^c  and  K2c  are  the  equilibrium  constants  for  the  formation 

of  the  carbamate,  the  first  and  second  dissociation  con¬ 
stants  for  H2C02  (carbonic  acid) ,  EA  ionization  constant, 
water  dissociation  constant,  the  inverted  Henry* s  constants 
for  CC>2  and  H2S  solubility  in  water,  the  first  and  second 
constants  for  H2S  dissociation.  The  average  coefficient 
for  the  ion  activity  is  given  the  symbol  y,  a  is  the  coef¬ 
ficient  for  the  water  activity  and  a  is  the  ratio  between 
the  un-ionized  ethanolamine  activity  and  the  activity  of 
water. 

Atwood  et  al.l  have  shown  that  for  the  system 
H2S-water-EA  the  form  of  the  dependence  of  the  average 
value  of  the  activity  coefficients  on  their  concentration 
can  be  taken  as  the  same  for  all  ions.  Klyamer  et  al.^ 
have  accepted  this  dependence  of  the  average  value  of 
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activity  coefficient  of  all  ions  on  the  ionic  concentration 
for  the  system  H2S-C02-water-EA  by  the  formulation  of  equa¬ 
tions  (5)  to  (8)  and  (10)  to  (12).  Atwood  et  al.^  obtained 
the  average  value  of  activity  coefficient,  y,  by  assuming 
that  the  amine  salts  are  similar  to  ammonium  salts.  They 
have  correlated  y  as  a  function  of  ionic  strength  of  the 
solution.  The  ionic  strength,  y,  is  defined  as: 


y  =  E  ^  m.  Z . ' 

li 


itu  =  molality  of  each  ion. 

=  electrical  charge  on  each  ion. 

The  values  for  the  dependence  of  the  correction  factor,  a,  on 

the  concentration  of  the  free  ethanolamine  and  the  values 

of  Kw,  Klc,  K2c,  K2y,  Hco  and  HR  g  in  the  range  of  temp- 

2  2 

eratures  25  to  120 °C  were  taken  from  the  work  of  Atwood 

et  al.^.  Values  of  K.  ,  K  and  K,.r  as  functions  of  temper- 

l  m  1Y  ^ 

11  12 

ature  have  been  taken  from  earlier  publications  ' 

The  analysis  of  equations  (1)  to  (13)  leads  to  the 
conclusion  that  the  most  important  parameters  which  enter 
into  many  reactions  are  the  so  called  degrees  of  sulphida¬ 
tion  (X  c)  and  carbonation  (X  ) .  It  has  been  shown 
H2b  C°2 

2-  + 

that  the  concentration  of  the  ions  S  and  H  can  be  neg- 

11  12 

lected  for  all  practical  cases  in  equations  (1)  to  (3) 

2-  - 

It  is  necessary  to  take  [CO^  ]  and  [OH  ]  into  account  only 


at  very  low  X  or  ^ 


co2  or  xh2s 


(i.e.  for  X 


,.r 


onpi 

r. 
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approximately  equal  to  0.01  mole/mole  of  EA) .  If  these 
simplifications  are  taken  into  account,  the  simultaneous 
solution  of  equations  (1)  to  (13)  allows  derivation  of 
expressions  for  the  calculation  of  the  equilibrium  partial 
pressures  of  C02  and  H2S  if  the  liquid  phase  composition 
and  temperature  are  known.  These  expressions  (as  derived 
in  Appendix  C)  are: 

K  Y2  A (A+B) 

P  =  - — -  •  —  •  -  - (14) 

2  Klc  Hh  g  a  a  (m-z-A-B) 


=  -L  .  y2  .  z (a+b) 

C0o  2  2  ,  ,  n.2 

2  K  a  a  (m-z-A-B) 
m  u 


(15) 


and 


A  =  m-z-B- 


K.  Kly  HCo: 


z 


a  (B-z) 


(16) 


where 


A 

B 


m  X 


m  X 


h2s  - 

Hh2s 

PH2s 

c02  ' 

Hco2 

PC02 

and 


z  =  [RNHCOO] 


Equations  (14)  to  (16)  have  three  unknowns;  PH  PCQ 

2 


and  z  ( [RNHCOO] ) .  Since  these  equations  are  implicit  in 
PH2S  and  Pco  [A  =  f(PH  g)  ;  B  =  f(PC0  )],  iterative 


fens 
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techniques  have  to  be  used  for  their  solution. 


For  a  degree  of  saturation  of  solutions  smaller 


than  0.7  moles/mole  of  EA,  the  physical  solubility  of 

and  CC>2  can  be  neglected  in  the  balance  equations^" '  ^ ^ . 

It  can  also  be  assumed  that  for  X  plus  X  ^  between 

H2b  C02 

0.05  and  0.7  moles/mole  of  EA,  CC>2  will  exist  as  RNHCOO 
and  HCO^  ions  and  as  HS  ions.  If  these  simplications 

are  taken  into  account,  equations  (14)  to  (16)  can  be  sim¬ 
plified  to: 

2 


K 


H2s  = 


W 


Ki  K1C  HH2S 


Y 
aa 


mXH2S  (XH2S  +  XC02} 


[hco31 


m 


+  1  xh2s  2xco2 


(17) 


Y 


C0o  v  22 

2  K  a  a 
m 


(XH2S  +  XC02' (XC02 


[hco3] 

m 


[HC03] 

m 


+  1  -  xH2s  -  2XC02 


(18) 


and 


[hco3]= 


m 


r  am  +  1  XH2S  2XC02 


K. 


KI  xco. 


am 


-=r  +  1  -  X„  c  -  2X 
am  H2S  CO2 


(19) 


where 


Ki  = 


Ki  kiy  hco. 


K  K 
w  m 
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Equations  (17)  and  (18)  can  be  further  simplified 

when  X  plus  X__  are  between  0.05  to  0.25  moles/mole 
H2S  C02 

of  EA  (concentrations  characteristic  of  the  regeneration 
process) .  For  such  cases  the  bicarbonate  ion  concentration 
can  be  neglected,  and  it  can  be  assumed  that  CO 2  exists 
only  as  carbamate  ion.  Therefore,  if  in  equations  (17)  and 
(18) ,  [HCO^]  =  0,  the  following  expressions  are  obtained: 

2 


K 


w 


”2S  Ki  K1C  HH2s 


mXH2S  (XH2S  +  XC02 
aa  (1  -  XR^S  -  2xcq2) 


(20) 


and 

P 


Y 


CO. 


„  2  2 

K  a  a 

m 


xco2  (xh2s  +  xco2) 

(1  ”  xh2s  “  2Xco2) 


(21) 


These  equations  are  much  simpler  than  equations  (14)  to 
(16)  or  (17)  and  (18)  and,  therefore,  can  be  used  for 
rapid  estimates  of  partial  pressures  of  the  acid  gases. 


If  only  H2S  or  C02  is  present  in  the  liquid  phase, 
equations  (14)  to  (16)  can  be  rearranged  to  give  an  expres¬ 
sion  for  Pj,  c  or  P--  .  If  X„  e  =  O,  i.e.  if  only  C0o  is 

11 2  <-(->2  ^2^  ™ 

present,  these  equations  reduce  to  the  following  form: 

1 


B  2K. 


'B  -  i  -  dg>  +  i~2  +  i '  B 


CO. 


K. 


A "  i(B  ■  i  ■ 


- (22) 


1  i  .  B  _  m 

2K2>  +  K,  2  B 


' 

■ 


i 
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where 


Ki  kiy  hco. 


K1  = 


K 


W 


aa 


Y 


and 


K2  = 


K  K 
m  w 


Ki  K1Y  HC0, 


a 


When  =0,  i.e.  if  only  is  present,  the  following 

2 

expression  for  calculation  of  P  is  obtained^. 

H2S 


m 


H2S  H  s(1-K") 


where 


t _ 

f  1-X„  _  ) 

/ 

H2s 

J 

2 

.  4 

+  K"  X 


h2s 


1-X 


h2s 


-  K"  X 


h2s 


-  (23) 


K 


K 


ii  _ 


w 


Kic  Ki 


Y 

aa 


Both  equations  (22)  and  (23)  are  implicit  and,  therefore, 
have  to  be  solved  by  iterative  techniques. 


Before  equations  (14)  to  (23)  can  be  solved  for  any 
temperature  and  any  loading  of  acid  gas(es)  all  the  equil¬ 
ibrium  constants  have  to  be  known  as  a  function  of  tempera¬ 
ture.  Also  the  correction  factor,  a,  as  a  function  of  free 
ethanolamine  concentration  and  average  ionic  activity  coef¬ 
ficient  (y)  as  a  function  of  ionic  strength  have  to  be 


I 
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1  11  12  13 

known.  In  the  literature  '  '  '  the  equilibrium  con¬ 

stants  as  functions  of  temperature;  the  correction  factor, 
a,  as  a  function  of  free  EA  concentration;  and  y  as  a  fun¬ 
ction  of  ionic  strength  have  been  presented  in  graphical 
form.  These  graphs  were  approximated  by  polynomials  using 
the  method  of  least  squares,  to  be  able  to  solve  equations 

(14)  to  (16)  by  using  a  digital  computer.  The  parameters 
for  the  polynomials  are  given  in  Appendix  H. 


The  Modified  Newton  Raphson  Method  was  used  to 

solve  the  above  equations.  For  each  iteration  a  better 

estimate  to  P  c,  P  and  z  was  obtained  by  the  following 

H2S  C02 

expressions : 


,i+l 

h2s 


=  P 


h2s 


>i+1 

ccu 


=  p 


CO, 


and 


i+1  i 
z  =  z 


af. 


ap 


h2s 


p1  p1 

h2s'  fco2' 


af. 


ap 


co. 


jpi+1  p1 

h2s'  co2' 


a  z 


pi+1  pi+1 
H?S'  C02' 


(a) 


(b) 


where  the  variables  superscripted  i+1  are  the  better  estim¬ 
ates  obtained  from  the  previous  ones,  superscripted  i. 
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Equations  (14)  to  (16)  are  rearranged  to  give: 

2 


K 


w 


fl  =  PH  C  - 

2  Ki  K1C  hh2s 


y  A (A+B) 


era  (m-z-A-B) 


f  =  P  _ L  .  _Y _  .  2  (A+B) 

2  C00  ^  22  ,  * 

2  K  a  a  (m-z-A-B) 
m 


f3  =  z  - 


a  K  K 
w  m 

Ki  kiy  hco. 


•  (B-z)  •  (m-z-A-B) 


The  partial  differentials  of  f^,  f 2  and  f^  with  respect  to 

Pu  c,  P^_  and  z,  respectively,  are: 

H2S  C02 


9  f . 


K 


9P 


=  1  + 


w 


H25 


Ki  Kic  hh2s 


(m-z-A-B) (zA+B)  +  (A+B)  H 


h2s 


(m-z-A-B) 


9  f , 


3p 


=  1  + 


Y 


ZH 


CO. 


£(m-z-A- 


B)  +  2 (m-z-A-B) (A+B) 


] 


CO. 


2  2 

K  a  a 

m 


(m-z-A-B) 


9  f . 


9  z 


=  1  + 


a  K  K 
w  m 

Ki  kxy  hco. 


(m  -  2z  -  A) 


Equations  (17)  to  (21)  were  used  to  obtain  the  initial  es¬ 


timates  of  P-.  0  and  P__  .  These  initial  values  along  with 

H2S  C02 


an  assumed  value  of  z  were  used  in  equations  (a)  to  (c)  to 
obtain  better  estimates.  Iteration  with  equations  (a)  to 
(c)  was  continued  until  the  difference  between  two  succes¬ 
sive  values  of  each  of  the  three  variables  was  less  than  1 


percent. 


' 


B.  Modifications  in  Klyamer  et  al.  Model: 


a)  Partial  pressures  were  substituted  by  fugacities 

in  equations  (9)  and  (13) .  This  was  considered  necessary 
especially  for  low  temperatures  and  high  total  pressures 
(pressures  greater  than  100  psia) .  The  fugacity  coefficient 
is  defined  as: 


K  = 


JK 

> 

K 


*pk  = 


fk 


where 


P^  =  partial  pressure  of  component  k 

f^.  =  fugacity  of  component  k 

(J)^  =  fugacity  coefficient  for  component  k. 

The  fugacity  coefficients  for  C02  and  H2S  were  calculated 
using  the  Modified  Redlich-Kwong  equation  of  state.  The 
expression  for  fugacity  coefficient  of  each  gas  (C02  or 

3 

H2S)  as  presented  by  Chueh  and  Prausnitz  is: 


m 


In  ({),  =  In 


v-b 


v-b 

ab. 


3/2  ,2 
RT  '  b 


yi  aik 

In 

v+b 

:3/2  b 

V 

v+b 

b 

Pv 

In - 

-  In 

— 

V 

v+b 

RT 

where 


. 
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a .  . 
ID 


2  2  *5 

(fi  .  +  ft  .)  R  T 
ai  a]  c .  . 

ii 


2  P 


c .  . 

ID 


c .  . 

ID 


z  R  T 
C  •  •  c  •  • 

_JJ _ ID 


V 


c .  . 
ID 


V  =  HIV1/3  +  V1/3) 
Ci  j  ci  cj 


c .  . 
ID 


(jO  •  +  0)  . 

0.291  -  0.08  (— - 1) 

2 


T  =  (T  T  )  2  (1  -  k  .  .  ) 

c .  .  c .  c  .  ij 

ID  ID 


The  molar  volume,  v,  is  that  of  the  gas  mixture,  and  is 
obtained  by  using  the  Redlich-Kwong  equation  of  state, 

p  =  RT  (a/T0,5) 

v”k  v(v+b) 

This  equation  is  rearranged  to  give: 

z3  -  z2  +  z (A-B-B2 )  -  AB  =  0 

where 
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and 


A  = 


aP 


2  2 
R  T 


B  = 


bP 

RT 


v 


zRT 

P 


The  cubic  equation  in  z  (the  compres: 


)  _L  JL  -L.  Y  -L_  / 


was  solved  and  the  largest  real  value  of  z  was  used  to  ob- 

2~  J2 


tain  v.  After  replacing  P  c  and  P„  by  the  respective 


fugacities  in  equations  (9) ,  (10)  and  (13) ,  the  set  of 

equations  (1)  to  (13)  was  solved  simultaneously  to  obtain 
the  fugacities  of  and  CC^.  These  fugacities  were  div¬ 

ided  by  the  fugacity  coefficients  of  each,  to  obtain  the 
respective  partial  pressures. 


b)  In  all  the  equations  used  to  calculate  partial 

pressure  of  acid  gases  (equations  (14)  ,  (15) ,  (17)  ,  (18)  , 

(20),  (21),  (22),  and  (23))  the  average  activity  coeffic- 

ient,  y,  appears  as  y  .  This  means  that  the  accuracy  in 

the  prediction  of  P  and  P  c  depends  on  the  accuracy  in 

C°2  H2b 

the  prediction  of  y.  As  mentioned  earlier  y  as  a  function 
of  ionic  strength  was  obtained  from  the  work  of  Atwood  et 
al.\  However,  when  the  solubility  data  for  the  system 
C02~H2S-water-EA  were  used  to  back-calculate  y  as  a  func¬ 
tion  of  ionic  strength,  the  simple  curves  reported  by  Atwood 
et  al.  were  not  obtained.  Instead  a  family  of  curves  was  ob¬ 
tained  depending  on  the  initial  EA  concentration  as  well  as 


■ 
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temperature.  Figure  8  shows  y  as  a  function  of  ionic 
strength  back-calculated  from  experimental  data  of  Lee, 

Otto  and  Mather^ '  ^ ^  and  also  of  Leibush  and  Shneerson^, 
for  different  temperatures  and  initial  DEA  concentrations. 
The  thick  line  shows  y  as  a  function  of  ionic  strength 
obtained  from  the  work  of  Atwood  et  al.\  As  is  evident 
from  this  figure,  y  is  not  only  a  function  of  ionic  strength, 
but  is  also  a  function  of  initial  EA  concentration  and  temp¬ 
erature.  It  can  also  be  observed  that  y  decreases  first 
and  then  increases  as  the  concentration  of  ions  in  the 
solution  increases  (i.e.  as  ionic  strength  increases) . 

With  increasing  initial  EA  concentration,  the  inflection 
point  moves  towards  higher  ionic  strength.  Theoretically 

this  type  of  behavior  is  predicted  by  the  Deby  e-Huokel 
5 

Theory  .  From  the  figure  it  can  also  be  observed  that  the 
temperature  dependence  of  y  for  the  same  initial  EA  concen¬ 
tration  is  small.  In  conclusion  it  can  be  said,  by  observ¬ 
ing  Figure  8  ,  that  the  simple  relationship  of  y  to  ionic 
strength  is  inadequate.  A  correlation  for  y  which  also 
incorporates  the  effects  of  temperature  and  initial  EA 
concentration,  is  therefore,  needed. 

In  this  work  an  empirical  function  of  temperature 
and  initial  EA  concentration  was  added  to  the  expression  for 
y  vs  ionic  strength.  The  expression  of  y  can  be  written  as: 


y  =  f  1  (y)  +  f2  (T  ,m) 


(XR  s  +  XCQ  )  >0.2  moles/ 

mole  of  EA 


. 

. 


■ 

. 


AVERAGE  ACTIVITY  COEFFICIENT 


IONIC  STRENGTH 

Figure  8  -  Activitv  Coefficient  vs  Ionic  Strength 
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Y  =  fx  (y) 


(Xh  S  +  XC0  ^  <0*2  moles/mole  of  EA 


where : 


y  =  ionic  strength 

T  =  temperature 

m  =  initial  EA  concentation . 


The  expressions  for  f^  and  are  given  in  Appendix  G. 

It  was  observed  that  the  prediction  of  experimental  data 
was  not  satisfactory  with  the  empirical  function,  when 

Xc0  plus  X^  s  was  less  than  0.2  moles/mole  of  EA.  The 
2  2 

reason  for  this  is  that  the  behavior  of  partial  pressures 

of  I^S  and  CO2  is  quite  different  at  these  low  loadings  of 

acid  gases.  Also  the  partial  pressure  of  f^S  does  not 

follow  the  same  trend  as  the  partial  pressure  of  C02  at 

these  conditions.  For  this  reason,  empirical  functions  of 

temperature  and  initial  EA  concentration  were  added  to  the 

expressions  for  predicting  P  *  and  P„  c  rather  than  to  the 

C02  H2S 

expression  for  y.  The  empirical  functions  are  presented  in 
Appendix  G. 

C.  The  Correlation  Based  on  Kent  and  Eisenberg^  Model: 


Kent  and  Eisenberg^  have  used  the  following  set  of 


chemical  equations  to  describe  the  I^S-water-EA  system: 


K 


eH+  +  RR'NH 


(1) 


(4) 


« 
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h2s 


H+  +  KS 


(6) 


HS 


+  S 


2- 


(7) 


h2s 


HH2S  ^H2s] 


(9) 


To  describe  the  C02~water-EA  system  the  equations  used  are: 


rr'nh2 

RR'NCOO 

h2o+co2 

h2o 

HCO" 


V 


H  +  RR'NH 


RR'NH  +  HCO" 


H  +  HCO. 


H+  +  OH' 


+  2- 
H  +  C03 


HC02[C02] 


(1) 

(2) 

(3) 

(4) 

(5) 

(8) 


RR'NH  represents  ethanolamine ,  where  R  stands  for  C2H^OH 
and  R'  for  H,  for  MEA;  and  R  and  R'  stand  for  C2H^OH  for 
DEA. 


The  second  set  of  equations  (for  C02-water-EA  system)  is 

4 

similar  to  that  proposed  by  Danckwerts  and  McNeil  .  Values 
of  equilibrium  constants  for  all  the  above  equations  have 
been  published*^.  However,  if  the  first  set  of  equations 
is  solved  for  H9S-water-EA  system  to  predict  P 


using  the 


. 
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published  equilibrium  constants,  an  unsatisfactory  fit  of 
the  published  vapor  data  results10.  Similarly  if  the  second 
set  of  equations  for  CC^-water-EA  system  is  solved  for  PCQ 

an  unsatisfactory  fit  of  the  experimental  data  results.  To 
solve  this  problem,  Kent  and  Eisenberg10,  treated  the  equil¬ 
ibrium  constant  ,  representing  the  EA  reaction  (for  the 
system  I^S-water-EA) ,  as  a  variable  and  accepted  all  the 
remaining  equilibrium  constants  as  published.  Using  pub¬ 
lished  vapor  pressure  data  for  15.3  weight  percent  MEA  sol¬ 
ution  and  20.5  weight  percent  DEA  solution,  those  values 
for  which  fit  the  data  were  found  for  each  temperature. 
For  CC^-water-EA  system  an  additional  equilibrium  constant, 

K 2,  was  used  because  C02  is  capable  of  reacting  directly 
with  either  MEA  or  DEA  to  form  the  amine  salt  of  a  substi¬ 
tuted  carbamic  acid  (the  carbamate  ion) .  To  solve  this  set 
of  equations,  only  K 2  was  allowed  to  vary.  For  the  ionic 
dissociation  reactions,  values  from  the  literature  were 
used.  For  ,  the  value  found  for  the  H^S-water-EA 

system  was  used.  Both  and  show  an  Arrhenius  depend- 

4.  10 
ency  on  temperature 

All  the  above  equations  were  combined  into  a  model 
for  the  mixed  system,  i.e.  the  system  containing  both  CC>2 
and  H2S  in  aqueous  EA  solution.  From  equations  (1)  to  (9) 
the  following  expressions  for  the  pseudo-equilibrium  con¬ 


stants  can  be  written: 


’ 


► 
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K1 : 

[H]  RR'NH 

-  (l) 

[RR'NH.  ] 

K2  1 

[RR'NH] [HCO~] 

=  - _A_  -  (ii) 

[RR'NCOO  ] 

K3  ' 

[H+] [HCO“] 

=  - -  (iii) 

[C02] 

K4  * 

=  [H+]  [0H_]  -  (iv) 

K5  * 

[H+] [CO,2- ] 

=  - r -  -  (v) 

[hco3] 

K6  ■ 

.  tH+HHS-]  -  (vi) 

[H2S] 

K7  = 

■  [H+]  [s2_]  -  (vii) 

[HS  ] 

hco2  = 

pco9 

=  ^  -  (viii) 

[C02] 

V 

ph9s 

=  — —  -  (ix) 

[H2SJ 

Also  the  following  balance  equations  can  be  written: 


m  = 

=  [RR'NH]  +  [RR'NH2+]  +  [RR'NCOO]  -  (x) 

“co2  ■ 

=  [R'RNCOO]  +  [C02]  +  [HC0~]  +  [CO^"]  (xi) 

mXH2S 

:  [HS~]  +[s2"]  +  [H2S]  - (xii ) 
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[RR'NH2+]  +  [H+]  =  [HC03]  +  [RR'NCOO]  +  2 [CO^  ] 

+  2  [S~ ]  +  [OH~ ]  +  [HS~] 


(xiii) 


where  all  the  concentrations  are  in  moles/liter  of  solution. 

Also  m  is  in  moles  of  EA/liter  of  solution.  Equations  (i) 

through  (xiii)  can  be  solved  simultaneously  to  obtain 

expressions  for  calculating  P_~  and  Pu  c  when  liquid  com- 

C02  H2b 

position  and  temperature  of  the  system  are  known.  These 
expressions  can  be  written  as: 


HH2s  A[H+]2 


H2S  K6K7  (1  +  Sp-) 

JA  « 


(14) 


H 


CO. 


CO. 


B [H+] 2 


K3K5  (1  + 


[H+]  m[H+]  . 

kc  kTkTk1”' 

5  2  5 


K 


A  (1  + 


B  (1  + 


(15) 


K2K5 


[H+] 


K?+  [H  ] 


— £  )  V  K2K5+K2[H+]  +  m[H+] 


■) 


K' 


(1  + 


m 


KXK 


r) 


(1  + 


m 


K.  K 


+ 


[H+]  (1  + 


(16) 


■ 
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where 


A  = 


XH2Sm 


h2s 


H 


h2s 


CO, 


B  = 


*C02- 


H 


CO, 


and 


+  PCO  K3 

K'  =  1  +  ■— — L  +  _ 

1  K2  HC0  ] 


Using  the  values  of  and  K 2  from  the  work  of  Kent  and 

Eisenberg^,  and  the  published  values  for  through  K^ , 

and  s,  equations  (14)  to  (16)  can  be  solved  for 
2  2 

PCQ  and  PH  s«  The  numerical  method  of  successive  approx¬ 


imations  was  found  to  be  adequate  to  solve  these  equations. 


D.  Comparison  of  Predictive  Capabilities  of  the  Models: 

The  starting  point,  for  both  Klyamer  et  al.  and 
Kent  and  Eisenberg  models,  is  the  postulation  of  chemical 
reactions  occuring  in  the  solution.  It  is  shown  in 
Appendix  D  that  for  the  system  H2S-C02-water-EA  only 
seven  independent  reactions  are  possible  and  that  the  sets 
of  reactions  postulated  in  the  two  models,  though  not 
identical  in  form,  are  chemically  equivalent.  The  model  of 
Klyamer  et  al.  is  more  sound  from  a  thermodynamic  stand¬ 
point  since  it  incorporates  the  effect  of  non-ideality  of 
the  liquid  phase  by  using  activity  coefficients.  The 
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substitution  of  partial  pressure  by  fugacity,  in  the  modi¬ 
fied  form  of  this  model,  lends  it  more  theoretical  support 
at  higher  pressures.  The  inclusion  of  fugacity  takes  care 
of  the  non-ideality  of  the  gas  phase.  The  model  of  Kent 
and  Eisenberg  incorporates  all  the  non-idealities  into  two 
parameters,  namely  and  This  changes  the  nature  of 

and  K 2  from  equilibrium  constants  to  two  empirical  par¬ 
ameters  which  enable  the  model  to  fit  the  experimental 
data  better. 


The  following  generalizations  can  be  made  about  the 
prediction  capabilities  of  the  three  models  (namely, 

Klyamer  et  al. ,  Kent  and  Eisenberg  and  Modified  Klyamer  et 
al.)  by  observing  Figures  9  to  26  in  Appendix  E  : 


(i)  In  general  for  MEA 
predicts  lower  values  of 

experimental  ones) ,  for  X 


,  the  model  of  Klyamer  et  al. 

P„  c  and  (compared  to  the 

H2b  C02 

plus  X  c  greater  than  ca. 

LUa  ri0o 


0.3  moles/mole  of  EA. 
mole  of  EA,  values  of  P 


For  loadings  lower  than  0.3  moles/ 


__  c  and  P__  higher  than  the  exper- 
H2S  C°2 


imental  ones  are  predicted.  The  prediction  worsens  with 
increasing  EA  concentration  and  acid  gas  loadings.  For 
DEA,  the  prediction  of  P  Q  and  Ppn  is,  in  general,  better 


than  that  for  MEA. 


.>  ' 
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(ii)  Only  for  simple  systems,  i.e.  systems  containing 

either  CC^  or  the  data  fit  obtained  by  Kent  and  Eisen- 

berg  model  is  fairly  good.  For  mixtures  (I^S-C^-water-EA 

system)  ,  on  the  other  hand,  prediction  of  P  and  P  c 

C02  H2S 

worsens  with  increasing  EA  concentration:  increasing  load¬ 
ings  of  CO2  and  H2S;  and  decreasing  temperature. 


(iii)  The  modified  Klyamer  et  al.  model  predicts  experi¬ 
mental  data  reasonably  well  over  most  conditions  of  temp¬ 
erature  and  EA  concentrations,  when  both  acid  gases  are 
present.  For  aqueous  MEA  solution  containing  either  CO ^  or 

H~S,  the  predicted  P  c  and  P^,-  values  are,  in  general, 

A  H2b  co2 

higher  than  the  experimental  ones.  For  DEA,  when  XCQ  or 

2 


X„  0  is  less  than  0.4  moles/mole  of  EA,  the  predicted  val- 
H2S 

ues  of  P^  and  P„  c  are  usually  slightly  higher  than  the 
C°2  H2b 

experimental  values.  However,  for  loadings  greater  than 

0.4  moles  (of  CO2  or  H2S)/mole  of  EA,  the  fit  is  reasonably 

good.  When  both  X^  c  and  X  are  less  than  0.1  mole/mole 

H2b  C°2 

of  EA,  this  model  reproduces  the  trend,  followed  by  exper¬ 


imental  P„  0  and  P~~  values,  however,  deviation  in  some 
H2S  c°2 

regions  is  still  present. 


By  comparing  the  predictions  of  each  of  the  three 
models  with  the  experimental  data,  the  following  conclusions 


can  be  drawn: 


■ 
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(1)  For  simple  systems  only,  i.e.  systems  containing 
either  C02  or  H^S,  the  model  of  Kent  and  Eisenberg  repro¬ 
duces  the  experimental  data  best,  over  most  conditions  of 
temperature,  EA  concentrations  and  acid  gas  loadings.  For 
the  same  systems,  in  most  cases  Klyamer  et  al.  model  also 
gives  a  reasonable  fit.  Modified  Klyamer  et  al .  model, 

on  the  other  hand,  usually  predicts  slightly  higher  PCQ 

2 

and  P  values,  compared  to  the  experimental  ones. 

H2b 

(2)  For  mixtures  of  C02  and  H2S  in  aqueous  solutions, 
the  Modified  Klyamer  et  al.  model  fits  the  experimental 
data  best,  especially  as  temperature  decreases;  EA  concen¬ 
tration  increases;  and  as  the  loading  of  acid  gases  in¬ 
creases.  For  loadings  of  CC>2  plus  H2S  greater  than  ca. 

0.8  moles/mole  of  EA,  the  modified  Klyamer  et  al .  model 

predicts  P  and  P„  c  which  are  slightly  lower  than  ex- 
C02  H2S 

perimental  values,  but  are  much  better  than  predictions 
from  the  other  two  models.  For  DEA  the  Kent  and  Eisenberg 
model,  though  not  better  than  modified  model  for  mixtures, 
is  better  than  the  Klyamer  et  al.  model  in  most  cases.  For 
MEA,  both  the  Klyamer  et  al.  and  the  Kent  and  Eisenberg, 
models  give  results  which  are  in  error  by  about  the  same 
order  of  magnitude.  Tables  I  and  II,  summarize  the  above 
observations  for  the  mixtures  of  C02  and  H2S  in  EA  solu¬ 
tions.  Each  entry  in  the  table  gives  the  percent  deviation 


in  the  predicted  partial  pressure  values  as  compared  to 
the  experimental  ones.  Models  I,  II,  and  III  stand  for 

Klyamer  et  al.,  Kent  and  Eisenberg  and  Modified  Klyamer 
et  al.  models,  respectively. 

(3)  Only  the  Modified  Klyamer  et  al.  model  is  able  to 

reproduce  the  trend  of  sharply  declining  Pu  c  and  P 

H2b  C02 

for  low  acid  gas  loadings  (when  both  and  X  Q  are 


less  than  0.1  moles/mole  of  EA) . 


,1.1  ,1  a  Is 

■ 


it  it  it 

TABLE  I:  PERCENT  DEVIATION  IN  PREDICTION  OF  PARTIAL  PRESSURES 
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SUMMARY  AND  CONCLUSIONS 


1.  An  experimental  apparatus  was  designed  and  con¬ 

structed  to  determine  very  low  partial  pressures  of 
CO2  and  at  high  temperatures. 

2.  The  apparatus  was  tested  by  reproducing  some  sol¬ 
ubility  data  of  H2S  and  C02  in  2.5  N  MEA  at  100°C.  A 
good  agreement  between  this  data  and  that  in  the  liter¬ 
ature,  confirmed  the  reliability  of  this  apparatus. 

3.  For  pure  CO2  in  2.5  N  MEA,  two  isotherms  at  80° 
and  100°C  were  extended  to  very  low  loadings  (ca.  0.02 
moles/mole  of  MEA).  Similarly,  for  pure  f^S  in  2.5  N  MEA 
the  80°  and  100°C  isotherms  were  extended  to  regions 

of  very  low  partial  pressures,  (about  0.0004  psia) 

4.  Measurements  were  also  made  for  mixtures  of  CC^  and 
H2S  in  5  N  MEA  at  100°C.  The  loading  of  each  CC^  and 
H2S  was  less  than  ca..  0.1  moles/mole  of  MEA.  No  pre¬ 
vious  data  were  available  for  this  region. 

5.  Of  the  two  theoretical  models  proposed  by  Klyamer 
et  al.  and  Kent  and  Eisenberg,  the  former  is  more  sound 
from  a  thermodynamic  standpoint  since  it  incorporates 
the  effect  of  non-ideality  of  the  liquid  phase  by  using 
activity  coefficients. 
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6.  The  modified  Klyamer  et  al.  model,  reproduces  the 
experimental  data  obtained  in  this  work,  as  well  as 
that  from  the  literature,  better  than  either  original 
Klyamer  et  al.  model  or  that  of  Kent  and  Eisenberg. 
(Figures  25  and  26) . 

7.  For  the  case  of  simple  systems,  i.e.  aqueous  -  EA 
systems  containing  either  CC>2  or  f^S,  the  model  of 
Kent  and  Eisenberg  gives  the  best  fit  of  experimental 
data. 
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APPENDIX  A 


CALCULATION  OF  P„  c,  P_^  ,  ETHANOLAMINE 

H2S  C02 

NORMALITY,  Xu  c  and  X__  . 

ri  ^  Jd  LUn 
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Calculations 


a)  Gas  Phase: 

Peak  area  for  each  gas,  N2'C°2  and  f^S  was  divided 
by  its  Relative  Molar  Response  (RMR)  factor.  Mole  percent 
of  each  acid  gas  was  obtained  by  the  formula: 

Mole  percent  H2S  = 

Peak  Area  H9S 
- —  x  100 

(RMR)H2s 

Peak  Area  H9S  Peak  Area  C09  Peak  Area  N9 

-  +  - -  +  - - 

(RMR)  (RMR)  (RMR).. 

H2b  C°2  N2 

Similarly, 


Mole  percent  C02 = 


Peak  Area  CO, 


(RMR) 


x  100 


H2S 


Peak  Area  H2S  +  Peak  Area  C02  +  Peak  Area  N2 


(RMR) 


h2s 


(RMR) 


CO, 


(RMR) 


N- 


RMR  factors  were  obtained  from  the  literature 


30 


RMR  for 


compound  i  is  defined  as 


Peak  area  of  i 


mole  percent  <J) 
x  -  x  100 


(RMR)  .  =  - 

mole  percent  of  i  Peak  area  of  (j) 
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where  <|)  refers  to  the  internal  standard  (benzene)  and  the 

factor  100  represents  the  RMR  of  benzene  arbitrarily  as- 

30 

signed  a  response  of  100  units  per  mole 


Since  the  total  pressure  was  always  near  atmospheric 

and  since  the  temperature  was  reasonably  high,  an  ideal  gas 
phase  was  assumed  in  calculating  the  partial  pressures  of 
E^S  and  C02.  Vapor  pressure  of  water  plus  amine  was  sub¬ 
tracted  from  the  absolute  pressure  in  B-3  to  obtain  the 
pressure  due  to  N2 ,  CO 2  and  H2S,  designated  as 

P (n2  +  h2s  +  co2) . 


Then: 


Partial  Pressure  of  H0S,  Pu  c 

Z  H2b 

Mole  percent  H2S 

=  - —  x  P(N9  +  co9  +  H?S) 

100  z  z  z 

Similarly, 

=  Partial  Pressure  of  C09 ,  P  n 

z  tu2 

Mole  percent  C09 

=  - -  x  P(N9  +  C09  +  H  S) 

100 


b)  Liquid  Phase: 

Normality  of  amine  (MEA) ,  NMEA: 

_  ml.  of  sulphuric  acid  used  x  normality  of  the  acid 

ml.  of  sample  used 


. 


Moles  CC>2  per  mole  of  amine,  XCQ  : 

ml.  of  hydrochloric  acid  x  normality  of  HC1 
used  to  titrate  BaCO^ 

2  x  N  x  ml.  of  sample  used 


Moles  of  H2S  per  mole  of  amine. 


(Normality  of  Iodine  x  \  ^  (ml.  of  thiosulphate 
.  initial  ml.  of  I0  used/  used)  x  normality  of 


MEA 


2  x  N 


x  ml.  of  sample  used 


APPENDIX  B 


EXPERIMENTAL  RAW  DATA 

a)  Experimental  Data  for  the  Solubility 
of  C02  in  MEA  Solutions. 

b)  Experimental  Data  for  the  Solubility 
of  H2S  in  MEA  Solutions. 

c)  Solubility  of  H2S-CC>2  Mixtures  in 
5  N  MEA  Solution  at  100°C. 
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(a)  Experimental  Data  for  the  Solubility  of  CC>2  in 


MEA  Solutions: 

(i)  2.5  N  MEA  Solution 


T/°C 


(Psia) 


XCQ  (moles  C02/mole  MEA) 


60 


80 


100 


0.0113 

0.012 

0.013 

0.0017 

0.00176 

0.0022 

0.003 

0.0027 

0.0097 

0.0149 

0.0126 

0.028 

0.033 

0.067 

0.071 

0.072 

0.081 

0.00017 

0.0002 

0.0005 

0.00053 

0.00063 

0.00067 

0.00076 

0.00098 

0.0010 

0.0052 

0.0135 

0.0141 

0.023 

0.109 


0.195 

0.194 

0.22 

0.031 

0.0313 

0.0325 

0.052 

0.052 

0.093 

0.112 

0.115 

0.158 

0.161 

0.214 

0.220 

0.217 

0.226 

0.0276 

0.0275 

0.029 

0.031 

0.0312 

0.028 

0.036 

0.035 

0.033 

0.0573 

0.080 

0.082 

0.094 

0.143 


'  ( 
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(b) 


(ii)  5.0  N  MEA  Solution 


T/°C 

PC02  (Psia) 

XC02  (Mole 

100 

0.00034 

0.035 

0.00053 

0.0375 

0.00077 

0.039 

0.0042 

0.060 

0.0081 

0.067 

0.016 

0.076 

0.043 

0.087 

0.054 

0.087 

0.130 

0.122 

0.190 

0.148 

Experimental 

Data  for  the 

Solubility  of 

Solutions : 

(i)  2.5  N  MEA  Solution 

T/°C 

PR  s  (Psia) 

XH  s  (Mole 
2 

80 

0.00483 

0.028 

0.0052 

0.027 

0.011 

0.041 

0.011 

0.042 

0.0166 

0.046 

0.0172 

0.0465 

0.023 

0.055 

0.025 

0.055 

0.05 

0.066 

0.065 

0.078 

0.079 

0.078 

0.081 

0.097 

0.085 

0.089 

0.162 

0.148 

CC^/mole  MEA) 


H2S  in  MEA 


H2 S/mole  MEA) 
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(c) 


T/°C 

PH  g  (Psia) 

Xu  c  (Mole  H0S/mole  MEA) 
H2b  1 

100 

0.0007 

0.012 

0.0012 

0.017 

0.0033 

0.017 

0.0081 

0.0216 

0.0164 

0.0294 

0.062 

0.053 

0.068 

0.048 

0.062 

0.053 

(ii) 

5.0  N  MEA  Solution: 

100 

0.00031 

0.026 

0.0051 

0.0344 

0.0097 

0.0303 

0.0156 

0.029 

0.040 

0.037 

0.055 

0.0435 

Solubility  of  H^S-CC^  Mixtures  in  5.0 

at  100°C.  (Partial  Pressures  in  Psia; 

acid  gas/mole  MEA) . 

P  p  y 

H2s  C02  H2s 

N  MEA  Solution 

X,  Mole  of 

XC02 

0.0019 

0.0015 

0.013 

0.0405 

0.0022 

0.0029 

0.0083 

0.0475 

0.0066 

0.0088 

0.010 

0.072 

0.026 

0.00031 

0.0264 

0.009 

0.047 

0.0024 

0.025 

0.026 

0.054 

0.0038 

0.0463 

0.023 

0.077 

0.071 

0.0133 

0.096 

0.085 

0.097 

0.013 

0.119 

0.097 

0.0021 

0.036 

0.012 

0.102 

0.048 

0.018 

0.078 

0.118 

0.164 

0.011 

0.116 

0.13 

0.063 

0.021 

0.082 

0.15 

0.0094 

0.032 

0.028 

0.16 

0.0013 

0.042 

0.006 

0.177 

0.088 

0.027 

0.083 

0.197 

0.056 

0.033 

0.06 

0.256 

0.00354 

0.053 

0.012 

. 

. 

> 


. 
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P  P  Y  Y 

H2S  C02  AH2S  C02 


0.305 

0.175 

0.022 

0.14 

0.44 

0.002 

0.06 

0.004 

0.46 

0.13 

0.035 

0.199 

0.47 

0.028 

0.053 

0.038 

0.490 

0.123 

0.046 

0.092 

0.557 

0.115 

0.037 

0.103 

0.60 

0.004 

0.084 

0.004 

0.61 

0.033 

0.066 

0.037 

0.64 

0.038 

0.063 

0.042 

0.64 

0.01 

0.08 

0.0124 

0.653 

0.015 

0.0796 

0.02 

0.716 

0.0132 

0.081 

0.0164 

0.76 

0.04 

0.069 

0.043 

0.93 

0.16 

0.071 

0.115 

0.988 

0.07 

0.071 

0.051 

1.04 

0.17 

0.069 

0.116 

1.  056 

0.069 

0.083 

0.057 

1.11 

0.006 

0.119 

0.006 

1.29 

0.0082 

0.128 

0.006 

1.35 

0.43 

0.064 

0.175 

* 

. 


APPENDIX  C 


DERIVATION  OF  EXPRESSIONS  FOR  CALCULATING 
P„  _  AND  P__  USING  THE  MODEL  OF  KLYAMER 

h2s  co2 

ET  AL. 
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Derivation  of  Expressions  for  Calculating  P  c  and  P  _ 

ii  2  ^ 


Using  the  Model  of  Klyamer  Et  A1 


[H2S] 

>H2S  = 


From  equation  7,  [H2S]  = 


Y2[HS  ]  [H+] 


K 


1C 


Therefore , 


p  =  Y2[HS  j  [H+] 

H2S  "  K1C  HH2S 

2-  11  12 

Neglecting  [S  ]  '  ,  from  equation  3 


[HS  ]  =  mXR  s  -  [H2S]  =  mXR  s  “  PH  S  HH  S  =  A 
2  2  2  2 


 y2  a[H+] 


H2s  K1C  HH2S 

aK 

From  equation  6,  [H  ]  =  - — — 

y2  [oh-] 


aK  A 
w 


H2S  K1CHH2S[0H'] 


(I) 


2-  +  11  12 

Neglecting  [S  ]  and  [H  ]  '  ,  equation  1  reduces  to: 


[RNH*]  =  [HC03]  +  [RNHCOO]  +  2 [C032  ]  +  [OH  ]  +  [HS  ] 


or 
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I 


[RNH3]  "  mXC02  -  PC02  HC02  +  A  +  !C032']  +  f0H'] 


=  B  +  A  +  [C032  ]  +  [OH  ] 


(i) 


substituting  [RNH^]  from  equation  5  into  equation  (i) 


K.aa  [RNH9] 

- - -  =  A  +  B  +  [CO  ]  +  [OH  ] 

Y  [OH  ]  J 


(ii) 


substituting  for  [RNH2 ]  from  equation  4  and  simplifying  for 


[RNH+] , 


K.aa  (m-z) 
[RNH^]  =  — 2-  -1-  ■_ - 


2r 


_ - j-  =  (A  +  B)  +  [CO"]  +  [OH  ] 

{y"  [OH  ]  +  K.aa  }  /'  / 


2- 

[CO^  ]  and  [OH  ]  are  considered  when  XCQ  is  less  than  0.01 
moles/mole  of  EA. 


From  the  above  equation; 


[OH  ]  = 


K^aa  (m-z-A-B) 
y2 (A  +  B) 


(iii) 


substituting  (iii) ,  into  equation  I 


K 


H2S  = 


W 


K1C  Ki  HH2s 


,  A (A  +  B) 
a  a  (m-z-A-B) 


From  equation  10, 


> 
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!  y2  [RNH+] 

p  zz  _  *  *  •  Z  •  _ 

CO-  K  2  2  rr).7tJ  ,2 

2  m  a  a  [RNH2J 

2 

,2  K.  aa 

1  Y  1 

—  -  •  1  -  •  2  •  --  -  -----  — 

Km  a2a2  y2 [OH-] [RNH2] 

-  (II) 

From  equation  4 , 

[RNH2]  =  m  -  z  -  [RNH*]  =  m-z-A-B  -  (iv) 

substituting  (iii)  and  (iv)  in  (II) : 

D  ^  _A_  •  ..r2  .  a  (a  +  B>  nn 

CO,  K  22  ,  2  (15) 

2  m  a  a  (m-z-A-B) 

From  equation  2, 

yl 

[ RNHCOO ]  =  z  =  mXc02  -  -  [HCO"]  -  ICo/'] 

=  B  -  [HCO~] 
or 

[HCO~  ]  =  (B  -  z)  -  (v) 

K  a[C09]  K1Y  HC09  PC09  [0H  ] 

[HCO_  ]  =  r-  —  =  - - - - - 

Y  [H  ]  K 

'  w 

substituting  value  of  [OH  ]  from  (iii) 

K1Y  HCO,  PCO,  K.  aa[RNH,] 

[HCO,]  =  - - - -  •  - -  -  (vi) 

K  Y  (A+B) 


'  ( 


. 
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From  equation  10, 


K 


m 


aV  [RNH2]2  P 
Y2  [RNH^ ] [RNHCOO] 


Multiplying  equation  (vi)  by 


a  [RNH2 ]  [RNHCOO] 


a  [RNH2 ]  [RNHCOO] 


(  =  1  ) 


and  rearranging, 


[rnh2]  = 


K1Y  HC02  Ki 


K  K 
w  m 


a  [HC03] 


kiy  hco2  Ki 

K  K 
w  m 


a  ( B- z ) 


( vii ) 


substituting  (vii)  into  (iv)  and  rearranging. 


Ki  kiy  hco 


A  =  m-z-B- 


2  . 


K  K 
w  m 


a  (B-z) 


(16) 


r* 


APPENDIX  D 


PROOF  THAT  CHEMICAL  EQUATIONS  POSTULATED 
IN  THE  TWO  THEORETICAL  MODELS  ARE  EQUIV¬ 
ALENT  . 
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PROOF  THAT  CHEMICAL  EQUATIONS  POSTULATED  IN  THE  TWO 


THEORETICAL  MODELS  ARE  EQUIVALENT 

32 

Zeggeren  and  Storey  have  shown  that  the  number 
of  independent  reactions  possible  for  a  given  system  is 
given  by: 

number  of  independent  reactions 

=  number  of  species  -  rank  of  atom  matrix. 

For  the  system  c02_H2^_water”EA' 
number  of  species 

=  12  (RR'NH,  RR'NCOO,  RR'NH*,  H20,  H+,  OH"" ,  H2S, 

HS_,  S2",  C02,  HCO~,  and  C032~) . 

Rank  of  atom  matrix  =  5  (N,C,H,0,S) . 

Therefore, 

number  of  independent  reactions  =  12  -  5  =  7 

Reaction  -  (1)  of  Kent  and  Eisenberg  model  is  a 
combination  of  reactions  (a)  and  (c)  of  Klyamer  et  al. 
model.  Similarly  reaction  (2)  of  Kent  and  Eisenberg  model 
is  a  combination  of  reactions  (a) ,  (c) ,  (e)  and  (f)  of 

Klyamer  et  al.  model.  Therefore,  though  the  seven  reactions 
postulated  in  the  two  models  are  different  in  form,  they 
are  chemically  equivalent. 
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APPENDIX  E 


COMPARISONS  OF  PREDICTED  PARTIAL  PRESSURES 
WITH  EXPERIMENTAL  DATA 


* 

-  Experimental  Data 

—  -  -  Prediction  by  Klyamer  et  al.  model 

—  •  -  Prediction  by  Kent  and  Eisenberg  model 

—  —  —  —  —  Prediction  by  Modified  Klyamer  et  al.  model 

of  Lee,  Otto  and  Mather  for  figures  9  to  24; 
of  this  work  for  figures  25  and  26. 
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Figure  9 


Comparisons  of  Predicted  Partial  Pressure 
of  for  2.5 

Published  Data. 


of  H2S  for  2.5  N  MEA  at  100°C,  with 
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Figure  10  -  Comparisons  of  Predicted  Partial  Pressure 

of  C02  for  2.5  N  MEA  at  100°C,  with 

Published  Data. 
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Published  Data. 
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(Psia) 
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Figure  12  -  Comparisons  of  Predicted  Partial  Pressure 

of  C02  for  2.5  N  MEA  at  40°C,  with 
Published  Data. 


(Psia) 


89 


Figure  13  -  Comparisons  of  Predicted  Partial  Pressure 

of  H2S  for  5  N  MEA  at  100°C,  with 

Published  Data. 


(Psia) 
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100.0  gt 


Figure  14  -  Comparisons  of  Predicted  Partial  Pressure 

of  C02  for  5  N  MEA  at  1Q0°C,  with 

Published  Data. 
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Figure  15  -  Comparisons  of  Predicted  Partial  Pressure 

of  H2S  for  5  N  MEA  at  40°C,  with  Published 

Data. 
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100.0 


Comparisons  of  Predicted  Partial  Pressure 
of  C02  for  5  N  M.EA  at  40°C,  with  Published 

Data. 


Figure  16 


(Psia) 
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Data. 


* 


(Psia) 
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Figure  18  -  Comparisons  of  Predicted  Partial  Pressure 

of  CC>2  for  2  N  DEA  at  100°C,  with 
Published  Data. 


(Psia) 
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Figure  19  -  Comparisons  of  Predicted  Partial  Pressure 

°f  for  2  N  DEA  at  50°C,  with  Published 

Data. 
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Figure  20  -  Comoarisons  of  Predicted  Partial  Pressure 

of  C02  for  2  N  DEA  at  50°C,  with 

Published  Data. 
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Figure  21  -  Comparisons  of  Predicted  Partial  Pressure 

of  H2S  for  3.5  N  DEA  at  100 °C,  with 
Published  Data. 
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of  C02  for  3.5  N  DEA  at  100°C,  with 
Published  Data. 
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Figure  23  -  Comparisons  of  Predicted  Partial  Pressure 

of  H2S  for  3.5  N  DEA  at  50°C,  with 

Published  Data. 


(Psia) 


100 


100.0 


10.0 


CN 

o 

u 

PM 


1.  0 


0.1 


0.01 


xco2 

Figure  24  -  Comparisons  of  Predicted  Partial  Pressures 

of  C02  for  3.5  N  DEA  at  50 °C,  with 

Published  Data. 
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Figure  25  -  Comparisons  of  Predicted  Partial  Pressure  of 

for  5  N  MEA  at  100°C,  with  Data  Obtained 


in  This  Study. 
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Pressure  of  CO^  for  5  N  MEA  at  100°Cr 
with  Data  Obtained  in  This  Study. 


APPENDIX  F 


COMPUTER  PROGRAMS 
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*****************  ******  ******  ******  *  ***  ******************** 
*******  *****************  *********************************** 


*  * 

*  PROGRAM  FOR  CALCULATION  OF  PARTIAL  PRESSURES  OF  H2S  AND  * 

*  C02  OVER  ET HANOLA  MINE  SOLUTIONS  — -  USING  MODIFIED  * 

*  KLYAMER  ET  AL.  MODEL  * 

*  * 


*  INPUT  VARIABLES  FOR  THE  PROGRAM  ARE  ,  * 

*  T  -  TEMPERATURE  IN  DEGREE  KELVIN  * 

*  MM  -  MOLES  OF  ETHANOLAMINE  PER  1000.0  GRAMS  OF  WATER  * 

*  KLMN  =  1  FOR  MONO-  AND  2  FOR  DI-ET HANOL AMINE  * 

*  * 


**************  *************************  ******************** 
*********************************************************** 


REAL  KIY,KI,KW,KIC,K2C,MA1,MA2,K1,KM, MU, MM,K2Y 
DIMENSION  DD  (3) 

DIMENSION  Y  (3)  ,PHI  (3) 

DO  101  1=1,3 
PHI  (I) =0.0 
101  CONTINUE 

READ  (5,1) T, MM, KLMN 

1  FORMAT  (2F15. 5,11) 

C  DETERMINATION  OF  VARIOUS  PARAMETERS  FOR  THE  ABOVE  READ  COND 

TT=1 000 . 0/T 

YY=0. 872 446E01- 0. 1 8331 2E01*TT+C. 33 752E00*TT**2 
KIY=EXP(-2.303*YY) 

YH=-0.450823E01+0.453045E-01*T-0.580136E-04*T**2 
HH2S  =  EXP  (-2. 303*YH) 

YC=-0. 248 352E01+0. 3611 75E-01*T-0. 4405 96E-04*T**2 
HC02=EXP (-2 . 303*YC) 

TW=  (T-273 .16) *1 . 8  +  32.0 

Y W=- 0. 75641 3+0. 27893 1E-0 1 *T W-G. 700 393E-04*TW **2+0. 747623E-0 
KW  =  EXP  (2.303*  (YW-15.0) ) 

YIC=0. 110835+0. 16671 8E-C1*TW-0. 6 84 879E- 04*TW** 2+0. 1069 72E-0 
KIC=EXP  (2.3  03*  (YIC-8.0) ) 

E2=-0. 29023 9E04/T+7 . 0880-0. 2379 E-01*T 
K2Y=EXP  (2 . 30  3*E2) 

K2C=0. IE-15 
IF  (KLMN-2)  2,3,3 

2  IF  (0.10E04/T-2.95)  4,7,7 

4  MA2=0.474  24E04*  (1 .0/T-0.  0025)  -3.6341 

K  M=EXP (MA2*2 . 303) 

GO  TO  10 

7  MA1  =  0.458E04*  (1 .0/T-0.  00295)  -1.50 

K M=EXP  (MA1*2.3C3) 

10  YIM=-2677 . 91/T-O . 61 69-0.4277E-03*T 

E3=EXP  (2. 303*YIM) 

KI=KW/E3 
AA=0. 8727 
GO  TO  20 

3  DA=0. 36  585E04*  (1 . O/T-0 . 0026)  -4.0 

KM=EXP  (DA*2 .303) 


■ 


. 


u  u 
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TTT=T-27  3 . 1 6 

YID=-1 .609+0.41 48 6E- 01* (TTT-37. 8) 

E33=EXP (YID) /0. 1E09 
KI  =  KW/E3  3 

AA  =  0. 857478+0.125767D-01 *MM+0. 4420 2D-C 2*MM** 2 

20  GAMA=0 . 50 
IF(KLMN-1)21,21,22 

21  WRITE(6,221) 

221  FORMAT  (1H1,/,1CX, •  ETHANOL  AMINE  PRESENT  IS  MEA*) 

GO  TO  18 

22  WRITE(6,222) 

222  FORMAT (1H1,/,10X,'  ETHANOL  AMINE  PRESENT  IS  DEA') 

18  WRITE (6,399) T, MM 

399  FORMAT (// ,2 X, *  TEMP.  OF  THE  SYSTEM= • ,F7. 2,//, 2X , • CONC. 

1 0F  ETHANOL  AMINE  IN  MOLES/1000.0  GRAM  OF  H20=',F5.2) 
WRITE  (6,412) HH2S,HC02, KIY,KIC,K2Y, K2C,KW,KM, KI 
412  FORMAT (//,2X, *  HENRY S  CONST.  FOR  H2S=* ,E1 4 . 5,//, 2X, ' HEN 
1RYS  CONST.  FOR  C02= ' ,El 4 . 5,///,2X , • EQIL.  CONSTANTS  K 
2IY,KIC,K2Y,K2C,KW,KM  AND  KI  RESPECTIVELY  ARE ' ,//, 2X, 7E 
316.6,////) 

DO  108  LLL=2 ,4,2 
DO  108  11=1,10,3 
XCO2=LLL/10 . 0 
XH2S=II/1 0. 0 

ALPHA=55.  51/(55. 51+MM*  (1 . 0  +  XCO2  +  XH2S)  ) 

INITIAL  VALUES  OF  PC02  AND  PH2  S  CAN  BE  OBTAINED  USING 

SIMPLIFIED  EQUATIONS 

IF  (1 .0-XH2S-2.0*XCO2)  68,  68,67 

67  K1=(KI*KIY*HC02)  /(KM*KW) 

CQ=K1/(AA*MM) +1. 0-XH2S-2. 0*XCO2 
QQ=XH2S+XC02 

HCM=SQRT  (  (CQ/2.0)  **2+ (KI *XC02) / (AA*MM) ) -CQ/2.0 
QT= (KW*GAMA**2*MM*XH2S*QQ)/(KI*KIC*HH2S*AA* ALPHA) 
PH2S=QT/(HCM+1. 0-XH2S-2. 0*XCO2) 

PC02=  (GAMA**2*XC02*QQ) / (KM* (AA*  ALPHA) **2* (HCM+ 1 . 0-XH2S 
1-2 . 0*XCO2)  **2) 

WRITE (6,415) PC02 , PH2  S 
GO  TO  220 

68  A1=MM*XH2S*0. 5 
B2=MM*XCO2*0 . 5 

IF  (XH2S) 416,416,417 

416  PC02=2. 0+B2/KM 
PH2S=0. 0 

GO  TO  418 

417  IF  (XC02)  419 , 419,420 

419  PH2S= (KW*A1**2)/(KIC*KI*HH2S*MM) * ( 1 . 0+XH2S/0 . 1 ) 

PCO2=0. 0 

GO  TO  418 

420  PH2S=(KW*A1* (A1+B2) ) /(KIC*KI*HH2S*MM) * (2.0+XH2S/0. 1) 
PCO2=2.0*  (Al+B2)/KM*  (1 .0+XCO2/0.1) 

418  WRITE  (6,415) PC02,PH2S 
415  FORMAT  (5X,E16. 5) 


■ 
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220  Z=0. 8*MM*XC02/(1 .0+4.0*XC02) 

231  ZZ=Z 

PHS=PH2S 
PCO=  PC02 
KL  =  0 

40  Al=MM*XH2S-HH2S*PH2S 
KL=KL+ 1 

B2=MM*XC02-HC02*PC02 
IF  (Al) 81,82,82 

82  IF  (B2) 81 , 83 ,83 

81  PH2S=0.95*PH2S 

PCO2=PCO2*0.95 
GO  TO  231 

83  FFK=MM-Z-A1-B2 
K  J=0 

IF  (FFK)  89,89,85 
89  FFK=MM/50 . 0 
85  HS=A 1 

12  RN=KI*AA*ALPHA**2 

OH=  (RN*FFK)/  (GAM A** 2*  (Al  +B2)  ) 

HPLUS= (KW*ALPHA) / (G AMA**2*OH) 

HC03=  (KIY*ALPHA*HC02*PC02) /  (GAM A**2*HPLUS) 

C03=  (K2Y*HC0  3)  /  (GAM A*HPLUS) 

RNH3=  (RN*  (MM-Z)  )  /  (G A MA** 2*OH+RN) 

S2=  (K2C*HS)  /  (GAMA*HPLUS) 

KJ=KJ+1 

IF  (KJ-10)  106,106,303 

106  MU  =  0 . 5* (Z+A1+HPLUS+OH+RNH3+4. 0*CO3 +HC03+4 . 0* S2 ) 

ETA=GAMA 

IF (MU) 300,300,301 

300  MU=2 . 0 

GO  TO  302 

301  IF (MU- 7.0) 302,302, 303 
303  GAMA=0 • 5 

GO  TO  93 

302  IF  (MU-0. 1 ) 74,74,75 

74  GAM A= 0.90230 8- 0.1995 77 E01*MU-0. 678868 E- 01*  (EXP  (MU) ) +0.05 
GO  TO  203 

75  IF  (KLMN-2)  54,55,55 

54  IF ( (XCO2+XH2S)-0. 2) 56, 56, 57 

56  GAMA=0. 399117-0. 338402 E- 02* MU+0 . 32 42 1 9E0 0*EXP (- MU) 

GO  TO  203 

57  GA  MA=0 .39  91  1 7-0 . 33 84 02E- 02*MU+0 . 32 421 9E00*EXP (- MU)  +0.1 
15*  (1 .0-  (0.21/(XCO2  +  XH2S)  )  **0.95)  +(  (XCO2  +  XH2S-0 . 4) /  (1 . 0 
2  +  EXP  (-  (7. 2-MM) )  )  )  *  (1 .1-EXP  (-  (374. 0-T)  /1 00. 0)  ) 

GO  TO  203 

55  IF  (  (XCO2+XH2S)-0. 2)  58,  58,59 

58  GAMA  =  0. 3664  0-0.  15975E- 01  * MU  +  0. 3443 79EC0* EXP (-MU) 

GO  TO  203 

59  IF  (XC02*XH2S) 60,61,60 

61  GAMA  =  0. 36640-0. 15975 E-01*MU+0.344379E0  0* EX P(-MU)  +  (0.10 
10*  (1 . 0- (0 .21/(XC02+XH2S)  )  **0.90)  +0.1*  (  (XCO2  +  XH2S-0 . 0  5) 


■ 
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2/  (XH2S+  XC02)  )*  (1 .05-EXP  {- (374. 0-T)  /I  00.0)  )  )*  (1.05-EXP( 
3-(5.22-MM)/10.0)  ) 

GO  TO  203 

60  GAMA=0. 3664 0-0. 15975 E- 01* MU+0.344379E00* EX P(-MU) +(0.10 

1  *  ( 1 . 0-  (0.21/(XCO2  +  XH2S) )  **0.90)  +0.5*  (  (XCO2  +  XH2S-0 .  30)/ 

2  (XH2S  +  XC02)  )  *  (1 .05-EXP  (-  (374 . 0-T)  /1 00 . 0)  )  ) 

203  IF (ABS ( (GAM A- ETA) /GAMA*1 00 . 0 ) -5 . 0)  93 , 93 , 1 2 

93  IF  (XC02) 91, 91,92 
92  IF  (XH2S)  94, 94,204 

91  PKK=  (KW*GAMA**2)/(KIC*KI*AA*ALPHA) 

PCC=SQRT ( ( (1 . 0-XH2S) /2.0) **2+PKK*XH2S) 

PHS=PH2S 

PH2S=MM/ (HH2S* (1 .0-PKK) ) * (PCC- (1 .0-XH2S) /2 . 0-PKK*XH2S) 

PCO=0.0 

GO  TO  933 

94  TK1= (KI*KIY*HC02*AA*ALPHA**2)/(KH*GAMA**2) 

TK2= (KM*KW*AA)/  (KI*KIY*HC02) 

B2=MM*XC02-PC02*HC02 

BB= (B2-MM/2. 0-0. 50/TK2) **2 
BM=MM/2. 0-B2 
PC  0=  PC  02 

PC02  =  B2/TK1*  (0.5  0/TK2-SQRT (BB+B2/TK2) +BM) / (0 . 5 0/TK2-SQ 
1RT  (BB  +  B2/TK2)  -BM) 

B2=MM*XC02-PC02*HC02 
PHS=0 . 0 
GO  TO  933 

204  FK 1 = (KW*GAKA**2) / (KIC*KI*HH2S*AA*ALPHA) 

FK2=  (GAMA**2)/(KM*  ( AA*ALPHA)  **2) 

FK3= (AA*KW*KM)/(KI*KIY*HC02) 

FFK=MM-Z-Al-B2 
IF (FFK) 953,953, 954 

953  FFK=MM/50 . 0 

954  Fl=PH2S-FKl*A1*  (A1  +  B2) /FFK 

NEW  VALUES  OF  PC02,PH2S  AND  Z  ARE  NOW  OBTAINED  USING 
MODIFIED  NEWTON- RAPHSON  METHOD 

DD  (1)  =1 . 0+FK1*  (FFK*  (2.  0*A1+B2)  +A1*  (A1+B2)  )  *HH2S/FFK**2 
PH2S  =  PHS- Fl/DD  (1) 

A1 =MM*XH2S-HH2S*PH2S 
FFK=MM-Z- A1-B2 
IF  (FFK)  955,955,956 

955  FFK=MM/50 . 0 

956  F2=PC02-FK2*Z* (A1+B2) /FFK**2 

DD (2)  =1 . 0+FK2*Z*HCO2* (FFK**2  +  2.0*FFK* (A1+B2) )/FFK**4 
PC02  =  PCO- F2/DD  (2) 

B2=MM*XC02-HC02*PC02 

FFK=MM-Z-A1-B2 

IF (FFK) 957,957, 958 

957  FFK=MM/50 . 0 

958  DD (3) =1. 0+FK3*  (FFK  +  B2-Z) 

525  F3=Z-FK3*  (B2-Z)  *FFK 

Z=ZZ-F3/DD(3) 

IF  (Z)  505,  505,506 


■ 
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505  ZZ=2.0*ZZ 
GO  TO  525 

506  IF  (  (XC02  +  XH2S)-0.2)  833,833,934 

934  DHL=ALOG  (ALPHA) /2.  303+8. 6  5-0.21  50E04/T 
PH20=EXP  (2. 303*DHL) 

P=PH2S+PC02+PH20 

PP=P/51 .70 

IF  (PP-50.0)  933,935,  935 

935  Y  (1) =PH2S/P 

Y  (2)  =  PC02/P 

Y  (3)  =PH20/P 

TRT=  (T-273. 16) *1.8  +  492.0 
CALL  FUGCF  (PP,TRT,Y,PHI) 

PH2S  =  PH2S/PHI  (1) 

PC 02= PC 02/P HI (2) 

GO  TO  933 

833  PPH=PH2S*  (1 . 0-EXP  (- XH2 S*  0 . 01 )  +2 . 0*EXP  (-  (0.06/XH2S)  )  +EX 
1P(-0.05/(XC02  +  XH2S*1 .0)  )  )  *(  1 . 0-0 . 0  0009/  (XH2S**2)  )  *  (1 . 0 
2+(XC02**2)/(2.5*XH2S)) 

PH2S=PPH 

PPC=PC02*  (1 . 0-EXP  (-XC02*0.d)  +2.0*EXP  (-  (0.07/XCO2)  )  +EX 
1 P (—0 . 05/(XCO2*0.2+XH2S) ) ) * (1 . 0-0 . 0 02/XH2S) *(1.0-0.004/ 
2XC02)  *  (1 . 0+XH2S**2/(1. 5*XC02)  ) 

PC02=PPC 

C  CHECKING  OF  CORRECTED  VALUES 

933  IF (ABS (  (PH2S-PHS) /PH2S*1 00.0) -5.0)  14,19,19 
14  IF (ABS ( ( PC02-PCO) /PC02*1  00. 0)-1.0)  16,19,19 
16  IF  (ABS  (  (Z-ZZ) /Z*100. 0)  -5 . 0)  69, 1  9, 1  9 
19  PHS=PH2S 
PCO=  PC02 
ZZ=Z 

IF  (KL-11)  40,40,405 

405  WRITE  (6,406) 

406  FORMAT (//,15X, 'NO  CONVERGENCE  AFTER  11  ITERATIONS'  //, 

1 1 0X , ' THE  LAST  ESTIMATES  WERE  AS  FOLLOWS  :',//) 

WRITE (6,1 00) XC02, XH2S, Z, PC02,PH2S 
GO  TO  108 

69  WRITE (6, 400) AA, ALPHA 

400  FORMAT (//, 2X , ' AA= ' , F7. 4,//,2X, *  ALPHA= ' , F7 . 4) 

WRITE  (6,77) KL 

77  FORMAT  (//,2X, 'NO.  OF  ITERATIONS  FOR  CONVERGENCE  TO  PAR 
1TIAL  PRESSURES=' ,13) 

PH2S=PH2S/51 .7 
PC02=PC02/51 .7 

WRITE (6,100) XC02 , XH2S, Z , PC02,PH2S 
100  FORMAT (//,9  X , ' XC0  2 ' , 13X, ' XH2S ' , 10X, 'RNHCOO (-)  '  ,6X, 

1 ' C02  PARTIAL  PRESS.  ' ,4X,  '  H2S  PARTIAL  PRESS .',//, 2X , 5E1 6 . 5 ,/ 
108  CONTINUE 
CALL  EXIT 
END 


. 


' 
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a*****************  **********************  **********  **  **  ****** 


*  * 

*  SUBROUTINE  FUGCF  CALCULATES  THE  FUGACITY  COEFFICIENTS  * 

*  OF  CARBON-DIOXIDE  AND  HYDROGEN-SULPHIDE  * 

*  * 

*  VALUES  IN  THE  CALLING  SEQUENCE  ARE  ,  * 

*  P  -  TOTAL  PRESSURE  IN  PSIA  * 

*  T  -  TEMPERATURE  IN  DEGREES  RANKINE 

*  Y  -  MOLE  FRACTIONS  OF  H2S,C02  AND  WATER  RESPECTIVELY  * 

*  PHI  -  INITIAL  FUG.  COEFF.  VECTOR  * 

*  * 


*********************************************************** 

SUBROUTINE  FUGCF  (P,T,Y,PHI) 

DIMENSION  Y  (3) 

DIMENSION  AIRKV  (3)  ,A  (4)  ,  Z  (3) 

DIMENSION  PHILN (3) , PHI (3) 

DIMENSION  ARKV  (3 , 3)  ,  BRKV  (3)  , PCIJV (3, 3)  ,TCIJV  (3,3) 

DIMENSION  GA  (3)  ,  GB  (3)  ,  PC  (3)  ,  TC  (3 )  ,  VC  (3) 

DATA  GA/0. 4340, 0.4470,  0. 421/ 

DATA  GB/0. 0882, 0.0911,0. 08368/ 

DATA  T C/6 72.7,240.0,1165.6/ 

DATA  VC/1.5606,2.3408,0.90/ 

DATA  PC/1  306.0,  1071.0,  3208.0/ 

PCIJV  (1 , 2 )  =0 . 57  41 E0  3 
PCIJV  (1,3)=0.2161 5E04 
PCIJV  (2,3)=0.10825E04 
TCIJV  ( 1 , 2) =0 . 369 7E0 3 
TCIJV  (1  ,3)  =0. 88  55E03 
TCIJV  (2,3) =0. 5289EC3 
R=1 0.73 
DO  200  1=1,3 

BRKV  (I)  =GB  (I)  *R*TC  (I)  /PC  (I) 

ARKV  (1,1) =G A  (I) *R**2 *TC (I) **2. 5/PC (I) 

DO  200  J= 1 , 3 

PCIJV  (J,I)  =PCIJV  (I,  J) 

TCIJV  (J,I)  =TCIJV  (I,  J) 

IF  (I-J)  2,200,2 

2  ARKV  (I,  J)  =  (  (GA  (I)  +GA  (J)  )  *R**2*TCIJV  (I,J)  **2. 5)  /(2.  0*  PCIJV  (I 
200  CONTINUE 

C  CALCULATION  OF  AMRKV, BRMKV 

AMRK V=0 . 0 
BM  RKV  =  0 . 0 
DO  120  1=1,3 
AIRKV (I) =0. 0 
BMRKV=BMRKV+Y  (I) *BRKV (I) 

DO  120  J=1 , 3 

AIRKV  (I)  =  AIRKV  (I)+Y  (J)  *ARKV  (I,  J) 

120  AMRKV=AMRKV  +  Y  (I)  *Y  (J)  *ARKV  (I,  J) 

C  CALCULATION  OF  VAPOR  MOLAR  VOLUME  OF  MIXTURE 

A  (1)  =1.0 
A  (2)  =-1.0 


■ 
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PBRT=P*BMRKV/  (R*T) 

ABRT= AMRKV/  (BMRKV* 1 0 . 7?*T**1 .5) 

A  (3)  =PBRT*  (ABRT-1  •  0-PBRT) 

A (4) =-ABRT*PBRT**2 
Z (1) =0. 0 
Z  (2)  =0.0 
Z  ( 3)  =0 . 0 
MTYPE=-1 

CALL  CUBEQ(A,MTYPE,Z) 

IF  (MTYPE)  130,140,140 
130  ZV=AMAX1  (Z  (1)  ,Z  (2)  ,Z  (3)  ) 

GO  TO  150 
140  ZV=Z  (1) 

150  V V=Z V*R*T/P 

C  FUGACITY  COEFFICIENTS  WITH  MODIFIED  R-K  EQN. 

QV VB=ALOG  (VV/ (VV-BMRKV)  ) 

Ql VB=1 . 0/  (VV-BMRKV) 

Q2RTB=2. 0/(10. 73 *T**1. 5* BMRKV) 

QV BV  =  ALOG  (  (VV+BMRKV) /VV) 

QARTB=AMRKV/  (10.73*T**1 . 5*BMRKV**2) 

QB VB=BMRK V/ (VV+BMRKV) 

DO  160  1=1,3 

PHILN (I) =QV VB+BRK V (I) *Q1 VB-AIRK V (I) *Q2RTB*QVBV+BRKV (I)  *QART 
1-QBVB)  -ALOG  (ZV) 

PHI(I)  =EXP(PHILN  (I)  ) 

160  CONTINUE 
RETURN 
END 


■ 
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*************************** ************ ********** ********** 
*  * 

*  SUBROUTINE  AMAXl  FINDS  THE  MAX.  Z  VALUE  * 

*  * 
*********************************************************** 

FUNCTION  AMAXl  (RR,SS,UU) 

AMAX1=RR 

IF (AMAXl- SS) 2,3 , 3 

2  AM AX1=SS 

3  IF  (AMAX1-UU)  4,5,5 

4  AMAXl =UU 

5  RETURN 
END 


1  *  >  * '  ‘  ■  *  •  *  1 
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*********************************************************** 


*  * 

*  THIS  SUBROUTINE  SOLVES  CUBIC  REDLICH-KWONG  EQUATION  * 

*  FOR  COMPRESSIBILITY  FACTORS  * 

*  * 


**  *  *  ***  **  **  **  *  ak*  *  *  ******  *  **  ^^nj:5{[^:^)(:^^^^.))i  **********  ********** 

SUBROUTINE  CUBEQ  ( A  ,  MTYPE, Z) 

DIMENSION  B  (3)  ,A  (4)  ,Z  (3) 

B(1)  =  A  (2)  /A  (1 ) 

B10V3  =  B  (1)  /3.0 
B  (2)  =  A  (3)  /A  (1) 

B  (3)  =  A  ( 4)  /A  (1 ) 

ALF  =  B  (2)  -  B  ( 1)  *B1 0V3 

BET  =  2.0*B10V3**3  -  B(2)*B10V3  ♦  B (3) 

BETOV  =  BET/2.0 
ALFOV  =  ALF/3.0 
CUAOV  =  ALFOV  **3 
SQBOV  =  BETOV  **2 
DEL  =  SQBOV  +  CUAOV 
IF  (DEL)  40,20,30 

20  MTYPE  =  0 

GAM  =  SQRT  (-ALFOV  ) 

IF  (BET)  22,22,21 

21  Z  ( 1 )  =  -2.0 *GAM  -B10V3 
Z  (2)  =  GAM  -B10V3 

Z  (3)  =  Z  (2) 

GO  TO  50 

22  Z (1)  =  2 . C*GAM  -B10V3 
Z  (2)  =  -GAM  -B1CV3 

Z  (3)  =  Z(2) 

GO  TO  50 

30  MTYPE  =  1 

EPS  =  SQRT  (DEL) 

TAU  =  -BETOV 
RCU=TAU+EPS 
SC  U=T AU-EPS 
SIR= 1 . 0 
SIS=1 . 0 

IF  (RCU)  31,32,32 

31  SI R=- 1.0 

32  IF  (SCU)  33,34,34 

33  SIS=-1.0 

34  R=SIR* (SIR*RCU) **0. 33333333 
S=SIS*  (SIS*SCU)  **0. 33333333 
Z  (1)  =  R  +  S  -  B10V3 

Z  (2)  =  -  (R  +  S)  /2 . 0  *  B10V3 
Z  (3)  =  0.86602540*  (R-S) 

GO  TO  50 
40  MTYPE  =  -1 

QUOT  =  SQBOV  /CUAOV 
ROOT  =  SQRT  (-QUOT) 


•  -  -  \  '  ■ 
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IF  (BET)  42,41,41 

41  PEI  =  (1.5707963  ♦  ATAN 
GO  TO  43 

42  PEI  =  ATAN  (SQRT  (1.0  - 

43  FACT  =  2 . 0 * S Q RT  (-ALF0  V 

Z ( 1 )  =  FACT*COS  (PEI)  - 

PEI=PEI+2. 0943951 

Z  (2)  =  FACT*COS  (PEI)  - 
PEI=PEI+2 .0943951 
Z  (3)  =  F ACT*COS  (PEI)  - 
50  CONTINUE 
RETURN 
END 


(ROOT  /  SQRT  (1.0 

ROOT**2 )  /  ROOT)  / 
) 

B10V3 

B10V3 

B10V3 


-  ROOT** 2)  )  )  /  3. 
'  3.0 
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APPENDIX  G 


EMPIRICAL  CORRECTION  FACTORS  FOR  MODIFIED 
KLYAMER  ET  AL.  MODEL. 
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CORRECTED  GAMA  FUNCTIONS  WHEN  XC02+XH2S  IS  >0.20 


(1)  FOR  MEA 

GAMA=0. 39911 7-0. 338402 E- 02* MU +0.32 42 19E00* EX P(- MU)  +0.1 
15*  (1.0-  (0.21/(XCO2  +  XH2S)  )  **0.95)  +  (  (XCO2  +  XH2S-0 . 4)  /  (1 . 0 
2  +  EXP  (-  (7. 2- MM)  )  )  )  *  (1 .1-EXP(-  (374. 0-T)  /1 00.0)  ) 

(2)  FOR  DEA 

(I)  WHEN  ONLY  ONE  ACID  GAS  IS  PRESENT 

GAMA=0. 3664 0-0. 15975 E- 01* MU+0. 3443 79E00*EXP(-MU)  +(0.10 
10*  (1 .0-  (0 . 2 1/  (XC02+XH2  S)  )  **0 . 90)  +0 . 1  *  ( (XCO2+XH2S-0 . 05) 
2/  ( XH2S  +  XC02 )  )  *  (1 .0  5-EXP  (  -  (374. 0-T)  /1 00.0)  )  )  *  (1 . 05- EXP  ( 
3- (5.22-MM) /10.0) ) 

(II)  FOR  MIXTURES  OF  H2S  AND  C02 

GAMA=0. 36640-0. 1597 5E- 01* MU+0. 3443 79E00*EXP(-MU) +(0.10 

1  * ( 1 . 0- (0.21/ (XC02  +  XH2S) )  **0.90) +0.5* ( (XC02+X H2S-0 . 30)  / 

2  (XH2S  +  XC02)  )  *  (1 .05- EXP  (-  (374. 0-T)  /1 00.0)  )  ) 


CORRECTED  PARTIAL  PRESS.  FUNCTIONS  FOR  XCO2+XH2S<0 . 20 


PPH  =  PH2S*  (1 . 0-EXP (-XH2S* 0.01 ) +2.0* EXP (- (0.06/XH2S)  )  +EX 
1P(-0.05/(XC02+XH2S*1 .0) ) ) * ( 1 . 0-0 . 0 0009/ (XH2S**2) )* (1.0 
2+ (XC02**2)/(2.5*XH2S)) 

PPC=PC02*  (1 . 0-EXP (-XCO2*0.01) +2.0* EXP (- (0.07/XC02)  )  +EX 
1P(-0.05/(XC02*0.2+XH2S))  )  *  (1 . 0-0 .0  02/XH2S)  *  (1 . 0-0. 004/ 
2XC02)  *  (1 . 0+XH2S**2/(1. 5*XC02)  ) 
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APPENDIX  H 


POLYNOMIAL  APPROXIMATIONS  OF 

Ki'  V  K1C'  K2C'  HH2S'  Km' 

K1Y'  K2Y'  HCO  '  Y  and  a 
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Kw  X  1()8 

ki  (DEA)  =  -1.609  +  0 . 0415 (T  -  37.8) 


T  =  °C 


For  MEA 


If  T  <  65°  K  =  e4-58  X  10Jd/T-0.00295)-1.5  .  ,K 

m 


If  T  >  65°C  K  =  e4’742  xlO  (l/T-0.0025)-3.634  T  =  °K 

m 


For  DEA 


v  3.66  x  10J (l/T-0.0026)-4.0  T  =  °K 

Km  =  e 


. 
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2  2 

Yn  =  a  +  a.x  +  a0x  +  a^x * 

1  o  1  2  3 

Y0  =  a^  +  a,x  +  a0e  x 

2  o  1  2 

Y3  =  ao/x  +  al  +  a2x 


Yi 

X 

a 

o 

al 

a2 

a3 

*1 

K1Y  =  e 

1/T  (K-1) 

8.724 

-1.833 

0.3375 

- 

HH2S  =  e  1 

T  (°K) 

-4.508 

0.0453 

-6  x  10"5 

- 

"co2  -  ^ 

T  ( °K) 

-2.4835 

0.03612 

4.4  x  10*5 

- 

«w  a  *Y1 

T  (°F) 

-0.7564 

0.0279 

7.0  x  10-5 

7.47  x 

IQ'8 

K1C  ■  21 

T  (°F) 

0.1108 

0.0167 

6.85  x  10"5 

1.07  x 

10~7 

a (DEA)  =  Yl 

m  (molality) 

0.8575 

1.258  x  10-2 

4.42  x  10"3 

- 

a (MEA)  =  Y 

m  3  0.4 

0.8727 

- 

- 

- 

Y  =  Y2 

M  ±  0.1 

0.9623 

-1.996 

-0.068 

- 

y(MEA)  =  y2 

p  >  0.1 

0.3991 

-3.384  x  10*3 

0.3242 

- 

y(DEA)  =  Y2 

p  >  0.1 

0.3664 

1.598  x  10‘2 

0.3444 

- 

Ki (MEA) 

Kw 

2.303  Y, 
e  3 

T  ( #K) 

-2677.9 

-0.6169 

4.277  x  10-4 

- 

